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Abstract 



We present the exact solution of a scalar field theory defined with noncommut- 
ing position and momentum variables. The model describes charged particles in a 
uniform magnetic field and with an interaction defined by the Groenewold-Moyal 
star-product. Explicit results are presented for all Green's functions in arbitrary 
even spacetime dimensionality. Various scaling limits of the field theory are analysed 
non-perturbatively and the renormalizability of each limit examined. A supersym- 
metric extension of the field theory is also constructed in which the supersymmetry 
transformations are parametrized by differential operators in an infinite-dimensional 
noncommutative algebra. 
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1 Introduction and Summary 



Despite the enormous amount of activity centered around the study of noncommuta- 
tive field theories (see P0-|3i for reviews), some of the most fundamental questions sur- 
rounding them have yet to be answered to full satisfaction. Foremost among these is 
whether or not these models can make sense as renormalizable, interacting quantum field 
theories. Non-local effects such as the interlacing of ultraviolet and infrared scales 
make conventional renormalization schemes, such as the Wilsonian approach, seemingly 
hopeless. In fact, this mixing is known to render some noncommutative field theories 
non- renormalizable, even if their commutative counterparts are renormalizable. A class 
of such models are asymptotically- free noncommutative field theories jSlEl; wherein the 
combination of asymptotic freedom and UV/IR mixing forces the couplings to run to 
their trivial Gaussian fixed points in the ultraviolet scaling limit. Such theories possess 
non-trivial interactions only when they are defined with a finite ultraviolet cutoff. This 
result is obtained by expansion about a translationally invariant vacuum of the quantum 
field theory. 

The analysis of UV/IR mixing and the renormalization of noncommutative scalar 
field theory to all orders of perturbation theory has been carried out in [ZI-[ni- For at 
least certain classes of noncommutative field theories, it may be that unusual effects such 
as UV/IR mixing are merely perturbative artifacts which disappear when resummed to 
all orders. This is of course a non-perturbative issue, which motivates the search for 
examples of exactly solvable noncommutative field theories. An example of such a model 
was presented and analysed in ilO|. It describes charged scalar fields on two-dimensional 
Euclidean spacetime in a constant external magnetic field and with a non-local quartic 
interaction. It may be thought of as the natural modification of a quantum field theory 
on a noncommutative space to one with noncommuting momentum space coordinates as 
well, i.e. as a quantum field theory on a noncommutative phase space. These models 
provide concrete examples in which one can see genuine quantum field theoretic issues, 
such as renormalization, at work in non-trivial exactly solvable theories. 

In jTU' the exact two-point Green's function of this model was obtained and the non- 
perturbative non-renormalizability of the field theory established. In this paper we will 
elaborate on this analysis and present various extensions and generalizations of this class 
of noncommutative field theories. In particular, we will extend the solution to arbitrary 
even dimensionality and higher order Green's functions, as well as presenting a detailed 
analysis of the various field theory limits which are possible in these models. In the course 
of this analysis, we will encounter new types of complex matrix models with external fields, 
for which we provide the necessary technical tools for analysis. Other aspects of these 
types of quantum field theories are discussed in [2]-^^, while their one-particle sectors 
are analysed in ^3^122] • Related exact calculations in nonrelativistic noncommutative 
field theory may be found in [23 121] ■ 
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1.1 Description of the Model 



The model we shall study in this paper is a particular deformation of scalar field theory 
in both position and momentum space. Consider complex scalar fields on even 

dimensional Euclidean spacetime M^" in an external constant magnetic field Bij, and 
with an interaction defined by replacing the usual pointwise multiplication of fields with 
the Groenewold-Moyal star-product, defined such that the components of the spatial 
coordinates x = {x^, . . . , x^") satisfy the relations 

x'i<x^ -x^ i<x' = 2ie'^ . (1.1) 

From (jl.lj) and the Baker- Campbell-Hausdorff formula it follows that the star-product of 
two plane waves is given by 

^ik-x ^ ^iq-x _ ^-ikid'^ Qj ^i(k+q)-x^ [1-2) 

Thus the interaction can be computed by expanding the fields in plane waves, i.e. by 
Fourier transformation, and (|1.2p defines the interaction part of the action for all Schwartz 
functions The background magnetic field, on the other hand, may be realized as 

the curvature of Landau momentum operators Pj = — i ^ — Bij x^ with 

[Pi,Pj] = -2iB,j . (1.3) 

The free part of the action is then defined by replacing the Laplacian operator of con- 
ventional scalar field theory with the Landau Hamiltonian P^. The 2n x 2n constant 
skew-symmetric matrices 6 = [6^^) and B = {Bij), characterizing respectively the non- 
commutativity of the coordinates and momenta, are assumed to be non-singular. 

It was shown in ^] that there is a natural regularization of these models which allows, 
in a simple manner, to make precise mathematical sense of the functional integrals defining 
the quantum field theory. Rather than expanding the boson fields in plane waves, it is 
more natural to expand them in the eigenfunctions of the Landau Hamiltonian P^ and 
thereby diagonalize the free part of the action. The Landau eigenfunctions 4>£^m{x) are 
labelled by two sets of quantum numbers £ = (£i, ...,£„), m = (mi,...,m„) G N", 
and the corresponding energy eigenvalues depend on only half these quantum numbers i. 
By expanding the fields as $(a;) = J2i m ^'rn,e. 4>e.,m{x) , the formal functional integration 
then becomes an integration over all complex expansion coefficients Mm,i- This suggests 
a regularization of the quantum field theory in which the collections of Landau quantum 
numbers (, and m are each restricted to a finite set which is in one-to-one correspondence 
with the set {1, 2, ... , iV}, where N < oo. With this regularization, the functional integral 
becomes a finite dimensional integral over N x N complex matrices. Similar matrix 
approximation schemes have also been studied in jHl^EHI^EH]- We remark that these 
matrix models of noncommutative field theory are not the same as those based on lattice 
regularization [^l3Uj. and in particular the scaling limits required in the two approaches 
are very different. 

As discussed in in the case of complex $^-theory deformed as described above, a 
finite matrix dimension N amounts to having both an ultraviolet and infrared cutoff at the 
same time, as there exist duality transformations which exchange infrared and ultraviolet 
divergences, while the action and all regularized Green's functions of the model are duality 
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invariant. In this paper we will carry this approach one step further and describe precisely 
how to remove the regularization, i.e. pass to the field theory limit N — >■ oo. In general, 
the interaction part of the action has a quite complicated form in the Landau basis (This is 
of course also true in the commutative case However, the quantum duality symmetry 
in this family of models suggests that they should be special when B = ±6~^ [lY. The 
crucial point is that at the special points B = ±6~^, the interactions acquire a simple 
form as well owing to a remarkable closure property of the Landau eigenfunctions under 
the star-product, and the field theory can be mapped exactly onto a matrix model |TH] . 
The significance of the points B = ±6~^ has been previously noted within different 
field theoretic contexts in [301^113 and within the context of noncommutative quantum 
mechanics in fT^ [W\ |^ . 

In this paper we will concentrate for the most part on extracting Green's functions of 
the noncommutative field theory in the limit where the regulator is removed. The latter 
limit corresponds to the large limit of the matrix model, and we can make it non-trivial 
simply by taking the 't Hooft limit of the matrix model |lOj. We thereby determine, in 
a completely non-perturbative way, how to renormalize the parameters of the model, i.e. 
give them a particular dependence on the regulator N such that the N ^ oo limit exists 
and is non-trivial. Even though this limit is essentially unique at the level of the matrix 
model, we will see that there is still some freedom left in defining the limit at the level of 
the field theoretic Green's functions. In ^U] the two-point Green's function was derived for 
one such limit in which it is translationally invariant. In the following we will use duality 
to argue that there is another limit, which mathematically is simply obtained from that 
of PU] by Fourier transformation, but which has a very different physical interpretation. 
We will use this result to argue that there may exist limits "in between" these other two 
scaling limits. 



1.2 Summary of Main Results 



From the mapping sketched above onto a finite-dimensional matrix model, we will derive a 
number of important consequences which lead to the exact solution of the noncommutative 
quantum field theory. Foremost among these features is the fact that the regularized 
field theory, for any interaction potential, is exactly solvable in the sense that there is a 
closed formula for the partition function at finite A^. In fact, we will establish the formal 
integrability of a more general matrix model which in the field theory amounts to having, 
in addition to the magnetic field, a confining electric field background for the charged 
scalars. In JH] this was established within the framework of the Hamiltonian quantization 
of the nonrelativistic, 2n-|- 1-dimensional fermion version of this field theory. It was shown 
that this model has a dynamical symmetry group GL{oo) x G'L(oo), with the free part of 
the model given as a linear superposition of Cartan elements of the representation and the 
interaction proportional to a tensor Casimir operator. The dynamical symmetry allows 
the explicit and straightforward construction of the eigenstates of the Hamiltonian, and 
some remnant of it should survive in the high-temperature limit of the corresponding 
partition function, which is essentially what we deal with here. In the following we will 
explore the integrability of this system from a different point of view. We will derive a 
determinant formula for the matrix integral and hence show that it is a tau-function of 
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the integrable, two-dimensional Toda lattice hierarchy. The GL(oo) x GL{oo) dynamical 
symmetry then presumably manifests itself as the rotational symmetry of the free fermion 
representation of the vertex operator construction of the Toda tau-function. The limit in 
which the electric potential is turned off reduces the partition function exactly to a KP 
tau-function. 

We will then study the large limit of the original matrix model. The Schwinger- 
Dyson equations of the matrix model succinctly sum up all divergences which one would 
encounter when solving the quantum field theory in perturbation theory. We shall com- 
pute the field theoretic Green's functions from the matrix model correlators, which also 
requires us to perform sums over Landau levels of certain combinations of the Landau 
eigenfunctions. The latter quantities capture the spacetime dependence of the Green's 
functions. As a warm-up and as a simple example where we can give complete results, we 
will first consider the limiting case in which the kinetic energy in the action is dropped. 
In this case we compute all Green's functions explicitly. 

We will then present the detailed computation of the two-point function G{x, y) = 
($^(a;) of the full model in various scaling limits of the quantum field theory with 

a complex <I>^ interaction potential. In the limit ^ oo with the quantity A^" = 
A^ ^/det{BJ4^^) held fixed, we find that the exact propagator admits the Fourier integral 
representation 

A_|_a-w.(a-(p^ + ^2)) e-(^-), (1.4) 

|p|<40FA 

where /i is the mass of the scalar field and W{X) is the propagator of the underlying 

matrix model which may be determined as the solution of the Schwinger-Dyson equations 
in the large A^ limit. The quantity A can be interpreted here as an ultraviolet cutoff, and 
in the limit A —>■ oo the function W reduces to 

_ ,1.5) 



with g the bare $^ coupling constant. This result was derived in |10j in the two- 
dimensional case. We will then argue that all higher Green's functions can be obtained 
from the two-point function by using Wick's theorem, i.e. all connected Green's func- 
tions in this particular scaling limit vanish, mainly as a consequence of the usual large 
A^ factorization in the matrix model. This type of scaling seems to be generic to the 
matrix regularization of noncommutative field theories [2^1 ED] ■ In particular, the contin- 
uum limits of the lattice derived matrix models of noncommutative field theory seem to 
be most naturally associated with the trivial Gaussian fixed point of the corresponding 
commutative theory [HJESl- Notice, however, that if we replace the bare coupling g by 
the renormalized coupling constant A^"~^5f, then the function A~^iy(A~^^) has a more 
complicated form reflecting the interacting nature of the quantum field theory, provided 
we simply interpret A as the finite mass scale set by the magnetic field. Then there is no 
trouble going to the full quantum field theory limit, and this suggests that an alternative 
physical interpretation for the infrared limit may also be given. This interpretation is 
consistent with the old folklore that spacetime noncommutativity provides a means of 
regulating ultraviolet divergences in quantum field theory. 
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We then exploit the duahty discussed above to show that there exists another scahng 
hmit of the quantum field theory, leading to a Green's function which is essentially the 
Fourier transform of ()1.4j) (up to rescalings). It corresponds to taking the limit N ^ oo 
with A^" = A^^^ a/ det(i?/47r) fixed and it produces an ultra-local propagator 

G^^{x,y) = -6'^^''\x-y)A~^w(^{A7rAf + A~^fi^'^, \x\<A-^/^, (1.6) 

where the parameter A now has the interpretation of an infrared cutoff. As remarked 
above, however, there is an alternative interpretation of both of these limits in which A 
may be interpreted simply as the mass scale set by the background magnetic field, i.e. 
A^" = a/ det(i?/47r). This follows from a remarkable scale invariance of the regularized 
Green's functions (and the underlying matrix model) which shows that, in the regularized 
quantum field theory, the length scale is set by the magnetic length /mag = det(-B)~^/'^"^. 
The two scaling limits thus described can be interpreted as zooming into short distances 
<C /mag in space in the infrared limit leading to ()1.4|) . and zooming out to large distances 
^ /mag in the ultraviolet limit leading to ()1.6|) . With this latter physical interpretation, 
one can proceed to analyse whether there exist well-defined intermediate limits with no 
zooming in or out of space at all. We thus find that the renormalization of the interacting 
noncommutative field theory depends crucially on the particular scaling limits that one 
takes, and also on their physical interpretations. 

We will also describe some other extensions of the present class of models which have 
the potential of providing exactly solvable, renormalizable, and interacting noncommuta- 
tive field theories. In particular, we present a detailed construction of a supersymmetric 
extension of the noncommutative scalar field theory in a background magnetic field. One 
of the most interesting aspects of these models that we shall find is that even at the level of 
the non-interacting theory, noncommutative geometry appears to play a role, in that the 
supersymmetry transformations are parametrized by elements of an infinite-dimensional, 
noncommutative algebra (or equivalently by infinite matrices). We use this observation to 
explicitly construct the most general naively renormalizable, noncommutative supersym- 
metric field theory, while it does not seem possible to construct local actions which possess 
this infinite parameter supersymmetry. This model thereby represents a supersymmetric 
extension of noncommutative field theory for which, like the scalar models studied in this 
paper, no commutative counterpart exists and noncommutativity is an intrinsic feature. 
We will not deal with the problem of actually solving any of these generalizations in this 
paper. Their exact solutions are left for future work which has the potential of deepening 
our insights into the basic structures underlying noncommutative quantum field theory. 

1.3 Outline 

In the ensuing sections we shall give detailed derivations of the results reported above, 
and describe many other interesting features of these noncommutative field theories. The 
structure of the remainder of this paper is as follows. In Section 2 we give the precise 
definition of the quantum field theory, describe its symmetries, and briefiy review the 
duality of In Section 3 we derive in detail, following 13 , the relationship between 
the quantum field theory and a large N matrix model, focusing our attention on the two- 
dimensional case for simplicity and ease of notation. We then use this mapping to describe 
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the nonperturbative regularization that will be employed for most of this paper. In 
Section 4 we prove that the regularized field theory is exactly solvable by deriving a closed 
formula for the finite N partition function, and relate it to the Toda lattice hierarchy. 
We also derive and explicitly solve the equations of motion in the large limit, and then 
evaluate the exact vacuum amplitude of the field theory. In Section 5 we study the Green's 
functions of the field theory in two dimensions, starting with those for which the free part 
of the action is dropped. All connected Green's functions are obtained in this limit and are 
invariant under Fourier transformation, as expected by duality. We then give a detailed 
derivation of the two-point function obtained in JU] in the infrared limit described above, 
and further argue that the connected correlation functions of the field theory are all 
trivial in this limit. In Section 6 we derive the generalizations of the preceeding results to 
arbitrary even spacetime dimensionality. In Section 7 the various other possible scaling 
limits of the quantum field theory are explored, including the infrared limit, and the 
different physical interpretations of the infrared and ultraviolet limits. We also briefly 
describe how non-trivial interactions can be obtained by (small) perturbations of the 
model. In Section 8 we present the detailed construction of the relativistic fermion version 
of the scalar noncommutative field theory, and combine it with the bosonic model to 
construct a supersymmetric, interacting noncommutative field theory in a background 
magnetic field. Appendix A describes an alternative interpretation of the model as a 
particular generalization of mean field theory for conventional (commutative) complex 
$^-theory, providing another motivation for the present work. In Appendix B we sketch 
a "stringy" interpretation of the duality of the model. The remaining Appendices C-G 
contain more technical details of the calculations. 



2 Formulation of the Model 



In this section we will start by defining precisely the noncommutative quantum field 
theory that we shall study in this paper. We shall also present a detailed analysis of the 
symmetries possessed by this model, and describe in what sense it may be regarded as an 
exactly solvable quantum field theory. 



2.1 Definitions 



We will consider the class of noncommutative scalar field theories defined by classical 
actions of the form 

S^ = j d^'^x $t(x) (-aD^-aD^ + zi^j $(a;) + V;($U$) (x) , (2.1) 
where a -|- 5" > and 



di - i Bij 
di + iBi^ x^ 



(2.2) 
(2.3) 



with di = d/dx^, i = l,...,2n. Here $(x) is a massive, complex scalar field of mass 
dimension n — 1 coupled minimally to a constant, background magnetic field proportional 
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to the 2n x 2n real, non-degenerate antisymmetric matrix B = (Bij), and in the back- 
ground of a confining electric potential proportional to (Bx)"^. Throughout we work on 
even- dimensional Euclidean spacetime M^"', and we use the notation = with 

the bar indicating complex conjugation. The interaction term is defined by a potential of 
the general form 

V{w) = J2fw\ (2.4) 



k>2 



and the subscript ^ means that the real scalar fields X[x) = are multiplied 

together in V^(X)(x) using the usual Groenewold-Moyal star-product, (X)^ = X -k X -k 
■ ■ ■ -k X {k times), where 

/ ^ f\x) = (27r)-^" j d'"A; j d'"g f{k) f\q) e e . (2.5) 

Here / is the Fourier transform of the field / with the convention 

Kk)= I d2"x e-'^--f{x) , (2.6) 



and for any two momentum space vectors k and q we have defined their skew-product by 

k-eq = e'^ kiqj . (2.7) 

We will assume that the noncommutativity parameter matrix 6 = (6^^) is non-degenerate. 
We also define 

x-By = B,j x' . (2.8) 



Most of the results we shall describe in the following will be obtained from the non- 
commutative complex $^ field theory with interaction potential given by 

V{w) = ^-w\ (2.9) 

and in the special instances when either a or a vanishes with a + a = 1. These cases 
correspond to charged particles in a background magnetic field alone. In Section I^TD we 
shall consider the case where both a = a = 0. The cases a,cr > will have the effect 
of removing the degeneracies of the Landau levels in the magnetic field problem. The 
special case a = 5" = | in ()2.H) corresponds to scalar fields in a harmonic oscillator po- 
tential alone and is closest to the conventional noncommutative field theories with no 
background magnetic field. The addition of such a harmonic oscillator potential to stan- 
dard noncommutative $^-theory is in fact necessary for its renormalization to all orders 
in perturbation theory j^. However, most of our analysis will focus on the instance where 
0" = 1, a = 0, and the interaction potential is given by (j2.9p . i.e. the noncommutative 
field theory with action 

S^ = j d^'^x <l>t(x) (-D2 + /i2)$(a;) + |$t^$^$t^$(3;) _ (2.10) 
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At the quantum level, the field theory is defined by the partition function 



D$ D<I>T e 



t ^-Sb 



:2.ii) 



and the non-vanishing Green's functions 



D$D<l>"f Y[<l>\xi)<l>{yj) e-^"^ , (2.12) 



1=1 



with the usual, formal functional integration measure D$ D$^^ = Yl^ d Re $(x) d Im 
The connected parts of the Green's functions ()2.12j) may be extracted from the generating 
functional which is defined by the functional integral with the coupling of the scalar fields 
to external sources J(x) and J{x), 



(2.13) 



We then have 



In 



Z^[J,J] 



fJ^5J{xi)5J{yi) 



:2.14) 



J=J=0 



To make the definition of the functional integral measure somewhat more precise we 
will expand the fields in a convenient basis of functions. Formally, one can choose any 
complete ortho normal set in the Hilbert space of L^-functions on R^". However, different 
choices of basis lead to different regularization prescriptions for the quantum field theory. 
For example, if one uses Fourier expansion of fields in a basis of plane waves e then 
the natural regularization imposed is the restriction of momenta k G M^" to the annulus 

< < A, where Aq is interpreted as an infrared cutoff and A as an ultraviolet cutoff. 
Removing the cutoffs then amounts to taking the limits Aq — > and A — »• cxd. This 
procedure is known to be rather difficult to carry out in noncommutative quantum field 
theory, and has in fact been done only at low orders in perturbation theory due to the 
complications set in by UV/IR mixing. In what follows we will circumvent such problems 
by using a different basis of functions that enables us to make sense of the quantum field 
theory non-perturbatively. 



2.2 Spacetime Symmetries 

Let us consider the spacetime symmetries of the field theory defined by the action ()2.H) . 
If n = 1, then the theory is rotationally invariant, because a constant magnetic field is 
rotationally invariant in two dimensions. However, parity is broken by the magnetic field. 
A parity transformation amounts to changing the sign of the magnetic field, B — ^ —B, or 
equivalently to interchanging a ^ a. What is much less obvious is that, in any dimension. 
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the model is also translationally invariant for a = 0. For this, we accompany a translation 
of the scalar field $(a;) by a gauge transformation in the constant magnetic field, 



^aix) = e ^{x + a) . (2.15) 

Note that the product on the right-hand side of ()2.15|) is just ordinary pointwise multi- 
plication of fields. Under ()2.15|) . the covariant derivative of $ transforms homogeneously, 

Di$„ = (Di$), . (2.16) 

An elementary calculation using ()2.5|) shows that the star-product also transforms in a 
simple way as 

$t^$^(x) = $U$(^x + (l-^5)a) . (2.17) 

It follows that the action (j2.H) with cr = is invariant under the substitution of $ by 
We will refer to this symmetry as "magnetic translation invariance" . 

If this symmetry persists at the quantum level, then the two-point Green's function 
G{x,y) = ($"''(x) $(?/)) assumes the form 

Gix,y)= e-'^'-^ygix-y) (2.18) 

with g{x) = G{x,0) a function only of one variable. This implies that the two-point 
function takes the same form in momentum space. Computing the Fourier transform of 
(jZTK|l gives 

G{p,q) = j d^^'x e'P-^ j d^"?/ e-'''-yG{x,y) 

= det(27rfi-i) e-'"'-^~'Pg(^B-\q-p)^ = det {2tiB-^) G {B'^q, B'^p) , 

(2.19) 

which for i? = 6'"^ is a special case of the duality to be discussed in Section 12.31 To 
interpret our later results it is important to note that the limits and B ^ oo are 

singular. At finite B the Green's functions in position and momentum space are both 
translationally invariant up to a phase. However, in the limit 5 — > only the translational 
symmetry in position space is maintained, with G{x,y) — > g{x — y). In momentum space 
this limit is ill-defined, but a direct computation yields G{p, q) = (27r)^" 5^'^"\p — q) g{p) 
and the momentum space Green's function becomes ultra-local in this limit. An analogous 
remark applies to the limit B ^ oo with the roles of G and G interchanged. 

One can now ask if the full symplectomorphism group of M^" is a symmetry of the field 
theory, as it is in the case of noncommutative gauge theories For this, we use the 
Weyl-Wigner correspondence to associate, for each Schwartz function f{x) on M?^ C, 
a corresponding compact Weyl operator / by 

(2.20) 
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where are the Hermitian generators of the Heisenberg algebra of n degrees of freedom 
which is defined by the commutation relations 



(2.21) 



We will represent this algebra on a separable Hilbert space Ti. With this correspondence, 
we can rewrite the action ()2.H) by using the usual rules 



= fg = f{x)g{x) , 
1 = /^ 

y"d2"x/(x) = Vdet(47r^) Tr^(/ 



(2.22) 
(2.23) 
(2.24) 

(2.25) 



where the dagger in ()2.23|1 denotes the Hilbert space adjoint and the trace of operators 
over Ti. in (j2.25|) is defined such that 



Tr 



n 



(e if-) = —2^1!^ 5(2") (p). 



By using (j2.22j) and (j2.24|) . along with the star-product identity 

* f{x) = x' fix) + i 6'^ djfix) , 
we may also derive the correspondence 



x'f = ^[x'f + fx' 



(2.26) 



(2.27) 



(2.28) 



It follows that the actions of the covariant derivatives in ()2.1|) on the Hilbert space H 
are given by 



55 = 



X-' 



(2.29) 



We see that the points where B 6 = ±11 have the special property that only one of the 
terms on the right-hand sides of the equalities in ()2.29|) is non-zero, and the thus the action 
assumes a particularly simple form. In particular, at the special point in parameter space 
where 



B9 = l , 

the action ()2.1|) is mapped under the above correspondence to 

= Vdet(47r0) Tr^ <!?\x)(^a (BxY + ij,'^^^{x) + a (!>\x)<!?{x){B 



[2.30) 



X] 



+ V(^<!>\x)<!>{ 



X 



(2.31) 
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and an analogous formula in the case B 6 = —H with the interchange a a. At these 
points, and only at these points, the action depends solely on the combinations 
and As we will see, this property makes the field theory exactly solvable at 

the full quantum level. Note that the representations of the two covariant derivatives in 
fl2.29|) also commute at these special points, so that the Hilbert space TC provides a Morita 
equivalence bimodule for the corresponding Heisenberg algebras that they generate. This 
equivalence is described from another point of view in Appendix B. 

The action ()2.31|) for a = possesses a large symmetry 

^{x) I — > U{x)^{x) , I — > ^\x)U{x)-\ (2.32) 

where U{x) is any invertible operator on Ti. In the spacetime picture, this GL{oo) sym- 
metry corresponds to the complex, fundamental star-gauge transformations 

$(x) I — > Ui<<^{x), <^\x) I — > <^{x)i.U*-\x) , (2.33) 

where U{x) is a star-invertible field on M^", 

U i.U'^-^x) = U*-^ i.U{x) = 1 . (2.34) 

With both a = a = the action has an even larger GL{oo) x GL(oo) symmetry 

$(x) I — > U{x)(!?{x)V{x)-\ ^\x) I — > V{x)^\x)U{x)-^ . (2.35) 

As we will see, this GL(oo) symmetry drastically simplifies the calculation of the Green's 
functions of the model. The unitary subgroup of this huge symmetry corresponds to 
invariance of the model under canonical transformations of M^" [36 . The symplecto- 
morphism W^i+oo-symmetry at the critical point ()2.30|) has also been observed within a 
different context in the corresponding noncommutative quantum mechanics pUl |22] • 

A somewhat different interpretation of the special point ()2.30p may also be given. The 
actions of the covariant momentum operators ()2.2j) and ()2.3|) may be written, by using 
(E2Z|), as 

Di^x) = + (dj<l>{x)+ iB'jj,x''i.<^{x)^ , 

b,^{x) = (l-BO)/ (^dj^{x)+ iB'^f^x''i<<^{x)^ , (2.36) 



where 



B' = i B , B" = i B . (2.37) 

1+B9 ' 1-B9 ^ ' 



This relationship is a very simple version of the Seiberg-Witten map between commutative 
and noncommutative descriptions of the same theory [3^. One can either work with the 
new background magnetic fields B' and B" using explicit noncommutative star-products 
in the definition of the Landau momentum operators, or else one can work with the 
original magnetic field B using the ordinary commutative product. The two descriptions 
are the same up to an overall rescaling. However, one of the transformations in (j2.37|) is 
always undefined at the points B = ±9~^, and in this sense they can be interpreted as 
points of "maximal" noncommutativity where the quantum field theory becomes exactly 
solvable.^ 



''Note that the transformations lfT!!7|l can be inverted to give B = B'{1 - 9 B'y^ and B = B"{\ 
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2.3 UV/IR Duality 



As demonstrated in |12j, the noncommutative field theory in the case of a quartic in- 
teraction also possesses a novel duality symmetry. As this duality will be important for 
the physical interpretations of the scaling limits of the regularized field theory that we 
will obtain, we shall now summarize the duality transformation rules, refering to ^2] for 
the technical details. The action ()2.10|) has a duality under Fourier transformation, i.e. 
it retains the same form when written in position or momentum space. Explicitly, it is 
invariant under the duality transformation 



det(5/27r) x) 



det(5i 



(2.38) 



with all other parameters of the model unchanged. This symmetry extends to the quantum 
level as a symmetry of rewriting the Green's functions (|2.12p in momentum space. 



G'b''(pi,...,Pr;gi,---,gr) = det (27r5 ^) 



(2.39) 



with xj = B^^pi, yj = B^^qj, I = l,...,r, and the transformations ()2.38|) implicitly 
understood on the right-hand side of ()2.39|) . The key to the proof of this fact is the 
existence of a regularization of the quantum field theory which respects the duality. This 
is essentially the regularization that will be exploited in the subsequent sections. Note 
that at the special points B = ±6~^, the field theory is essentially "self-dual" in that it 
is completely invariant under Fourier transformation, and the duality identifies the two 
special points. This is consistent with what was found in the previous subsection. In 
Appendix B we point out an alternative "stringy" interpretation of this duality at the 
special point. 



3 Matrix Model Representation 

We are interested in expressing the action ()2.1|) with respect to the expansion of functions 
on M^" in a suitable basis. To simplify the presentation of our results we will first focus 
on the two-dimensional case n = 1. The generalizations of our results to arbitrary even 
dimensionality is remarkably simple and will be described later on in Section IHl In this 
section, and until Section IHl we define 9 = 9^^ and B = B12, both of which are assumed 
to be positive parameters. We will begin by deriving the matrix form of the action (j2.H) 
for n = 1, and then use this representation to define the formal functional integral of 
the quantum field theory as an integral over infinite-dimensional matrices. This will then 
enable us to describe a natural non-perturbative regularization of the model. 

9 B")~^, in which the roles of the two special points B' — +9^^ and B" = —9^^ are interchanged. The 
Seiberg-Witten transformation then asserts that the two points are equivalent, in that they lead to the 
same noncommutative quantum field theory. 
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3.1 Mapping to a Matrix Model 



Let us begin by introducing the complex coordinates 



Z = X + IX , Z = X —IX 



(3.1) 



and the ladder operators 

1 / 



2 V Ve 
2 V Ve 



Ve 



bt = ^ (-V9d + ^ 

2 V Ve 



(3.2) 



Here d = di— 182 and d = di+ 182- The operators ()H.2|1 generate two commuting copies 
of the usual harmonic oscillator commutation relations 



[a, at] = [b,bt] = 1 
[a, b] = [a, bt] = 



(3.3) 



The Landau states are defined by 



\i,m) 



(bt)"^-i 



I vac) 



(3.4) 



where |vac) is the two-particle Fock vacuum with a|vac) = b|vac) = 0, and £,m are 
positive integers. The normalized Landau wavefunctions in position space are given by 



4>e,m{x) = {x\i,m) . 



(3.5) 



The ground state wavefunction (f)o{x) = (pi,i{x) can be computed explicitly by solving the 
differential equations a 0o = and b 0o = using (j3.2p to get 



-\z\^/2e 



(3.6) 



'TTt 



while in general they can be computed by introducing the generating function 

-1 ^m~l 



VsAx) = J2 



4>e,m{x) 



1 VW^-Vi^^^- 

Using ()3.2|l - ()3.f)j) . an elementary calculation gives (see for example jTHj, Appendix A) 



(3.7) 



'ixe 



and from ()3.8|) the Landau wavefunctions may be extracted as 



1 



1 am— 1 



v/(£-l)!(m-l)! ds^~^ df 
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-^sAx] 



(3.9) 



s=t=0 



The advantage of working with the functions (pe^mix) is two-fold. First of all, using 
the generating function ()3.8p it is easy to check that they form a complete orthonormal 
basis of one-particle wavefunctions, and 



d X (pi^rn * (j)i',m'{x) = 5ei' 5mm' (3.10) 



with (pi^rn = <Pm/- Secoudly, a straightforward computation using ()3.8|) and the represen- 
tation ()2.5p yields the identity (see for example jl3j, Appendix A) 

Vs,t^Vs',t' = e'''Vs,t' , (3.11) 

which is equivalent to the star-product projector relation 

4'e,m'*' 4'e',m' = ^j-^-Q ■ (3-12) 

This relation of course just reflects the well-known fact that the functions VAttO (j)e^rn 
provide the Wigner representation of the rank 1 Fock space operators (pi^m = 

These facts suggest an expansion of the scalar fields in the action ()2.H) with respect 
to the Landau basis as 

oo oo 

<i>(x) = J2 ^(^) ' '^^(^) = Yl ^ '^^'™(^) • (3.13) 

£,m=l t,m=l 

It is natural to assemble the dimensionless complex numbers into an infinite, square 
complex matrix M = {Mim)i,m>i- To ensure convergence of the expansion ()3.13|) . we 
will regard M as a compact operator acting on the separable Hilbert space Ti. = S{N) 
of Schwartz sequences {am)m>i with sufficiently rapid decrease as m — > oo. Then, the 
projector relation ()3.12|1 implies that the star-product of fields $(x) and $'(x) corresponds 
to ordinary matrix multiplication [25j, 

oo 

e,m=i 

while ()3.12|) along with the orthogonality relation ()3.10|) implies that spacetime averages 
of fields correspond to traces over the Hilbert space 7i, 

/oo 
d^x $(x) = Aire J2 ^n^m = 47r0 Tr^(M) . (3.15) 
m=l 

The star-product powers ($^ interaction potential V^($^ 'k^){x) may thereby 

be succinctly written as traces of the matrices (M'^M)''. For the kinetic energy terms in 
()2.1|) . we use the definitions ()2.2|) and ()3.2p to write 



L + 50)^(ata + 0+(l-i?0)^ 
+ {B^e^-l) (atbt + ab) 



(3.16) 
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along with the usual harmonic oscillator creation and annihilation relations applied to the 
Landau wavefunctions ()3.5|) to get 

a (pi^rn = - 1 4>e-l,m , 3^ 4>e,m = 0£+l,m , 

h(j)e^rn = Vm - 1 , b"^ 4>e,m = y/m ■ (3-17) 

The operator is obtained from ()3.1(jj) by a reflection of the magnetic field B —>■ —B. 

We thereby find that the total action ()2.H) may be expressed as the infinite-dimensional 
complex matrix model action 

5b = {a + d)[B'9^-i) Tr^(rLMtr^M + MtrLMrJ 

+ [a + d + {a - d)B'^e'^^ Tr^ (M£Mt) 
+ (^a + a - (a - a)B^e^^ Tr ^ {M^ £ M) 

+ ATTen^ Tr^ (M^M) + A-nO Ti [M^ M) , (3.18) 
where we have introduced the diagonal matrix £ = {£f,rn) with 

£i^ = Ati(i-^ 5,m (3.19) 



2^ 

and the infinite shift matrix T^o = (Xoo/m) with 

v/47r(m - 1) Se,m-i • (3.20) 



oo,£m 



Note that the noncommutativity parameter ^ now only appears as an explicit coupling 
parameter in the action ()3.18|1 . The noncommutativity of spacetime has simply become 
the noncommutativity of matrix multiplication. The first three traces in the action ()3.18j) 
are similar to the clock and hopping terms that appear as the kinetic parts in the lattice 
derived matrix models of noncommutative scalar field theory PU] . 

The maximal symmetry group of area-preserving diffeomorphisms, which appears at 
the special points in parameter space where either cr = or a = 0, and B6 = ±1, is 
particularly transparent within this formalism. At these points, the first term in ()3.18j) 
vanishes and the resulting action depends only on the combinations M'^'M and MM^". In 
particular, for 5" = it possesses the GL(oo) symmetry 

Ml — > U -M , M"^ I — > -U-^ . (3.21) 

As we will discuss below, in this case the noncommutative quantum field theory has the 
same fixed point as the complex one-matrix model in the 't Hooft limit. Note also that 
in this case the Landau states diagonalize the Landau Hamiltonian operator D^, with the 
standard harmonic oscillator spectrum 

-D20,,„ = 45(^£-0 (3.22) 

The GL(oo) symmetry can thereby be physically understood as a consequence of the in- 
finite degeneracies of the Landau levels ^, and it acts through rotations of the magnetic 
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quantum numbers m. The phase space now becomes degenerate and the Landau wave- 
functions depend on only half of the position space coordinates, leading to a reduction of 
the quantum Hilbert space a.t B = 9~^. Henceforth we shall deal only with the quantum 
field theory defined at this special point. Note that from ()3.16|) it follows that while the 
operator at B9 = 1 is given by the a,a^ oscillators, the Hamiltonian is given in 
terms of the commutant b, oscillators and the eigenvalue equation ()3.22|) is modified to 



3.2 Regularization 

Let us consider now the quantum field theory in the matrix model representation ()3.18j) . 
whereby the scalar field variables can be identified with infinite matrices. The Jacobian 
of this transformation is formally 1, and the partition function (j2.1H) is then given by the 
functional integral for the N = oo complex matrix model. With this rewriting we can 
now properly define the functional integration measure appearing in 1)2.111) . Namely, we 
replace the infinite matrix variables by finite N x N matrices, and then take the large 
N limit of the matrix model. In other words, we restrict the quantum numbers of the 
Landau wavefunctions to £, m = 1, . . . , with N < oo. There are several advantages 
to this matrix regularization. For example, as discussed in [12 , a finite matrix rank 
provides both a short distance and low momentum cutoff simultaneously, and thereby 
avoids the mixing of ultraviolet and infrared divergences which seem to make the analysis 
of noncommutative perturbation theory to all orders hopeless. Introducing an ultraviolet 
cutoff alone in the matrix form is not enough to regulate all possible divergences. Ul- 
traviolet and infrared divergences cannot be clearly separated and one needs to regulate 
them both at the same time. This point will be further elucidated in Section [7| In this 
respect, the Landau basis for the expansion of noncommutative fields is far superior to 
the conventional Fourier basis, in that it naturally leads to a non-perturbative regular- 
ization which avoids dealing with the usual difficulties of noncommutative quantum field 
theory. It is important to note though that this matrix regularization does not impose 
any a priori restrictions on the allowed noncommutativity parameters 6, in contrast to 
what occurs in lattice regularization |3Uj . 

However, while the matrix regularization is close to a cutoff in the kinetic energy term 
of the action ()2.1|) . cannot itself be an ultraviolet cutoff, because it is not determined by 
any mass scales of the theory. Without a dimensionful scale there can be no dimensional 
transmutation of the coupling constants of the model, which is the essence of renormaliza- 
tion in quantum field theory. The renormalization ensures finiteness of physical quantities 
and requires some set of observables to define renormalization conditions. The way to 
define the theory is explained in [21] and consists of first introducing the kinetic energy 
cutoff, then taking the limit 6 oo with the cutoff scaled such that it remains finite, 
and finally removing the cutoff. This limit is equivalent to the conventional large limit 
of the complex matrix model |2Zl EHI- In other words, we should take the limit N ^ oo 
while keeping fixed the dimensionful ratio 




(3.23) 



N 



(3.24) 
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which simply defines the large 9 limit of the model ()2.H) . The true ultraviolet cutoff ()3.24|) 
has a very natural physical interpretation. At -B = 9~^, it is associated with the energy 
2B{2N — 1) = IGttA^ — 2B of the A^**^ Landau level, which remains finite in the large A^ 
limit just described. Since B —>■ a.s N —>■ oo, the spacing between Landau levels also 
vanishes. Thus taking the limit described above is equivalent to filling the finite energy 
interval [0, IGttA^] with infinitely many Landau levels and an infinite density of states. 
This limit can be thereby regarded as a combination of those found in the corresponding 
noncommutative quantum mechanics in jl5|, 1^ 1^ . whereby the Laughlin theory of the 
lowest Landau level is recovered, and in ^ni, where all Landau levels survive. 

The large 6 limit of the model ()2.1|) thereby defines a field theory with a finite cutoff. 
The quantum field theory in this limit coincides with the 't Hooft limit of the matrix 
model and, as we will see below, is exactly solvable (at the special point where B6 = 1). 
At this stage, we should renormalize the regulated quantum field theory, i.e. fine-tune 
the parameters of the model in such a way that the limit A ^ oo is well-defined and non- 
trivial. For this, we rescale the massive coupling constants /x^ and gk to the dimensionless 
ones 



2 

= ^ , 9k = j^. (3.25) 



Let us now specialize to the case a = 1, a = in (|2.im) . The functional integral (j2.1H) 
may then be properly defined through the large A^ limit of the matrix integral 

Zi,{E) = yDMDMt e-^T'^(*^^*^'+^(^^'^^)) , 

DM DMt ^ fr dM,^ dM,^ ^ ^^ ^6) 

£,m=l 

with the integration extending over the finite-dimensional linear space oi N x N complex 
matrices M = {M(,m)i<£,m<N, and 

n^) = E|^' (3.27) 



k>2 



the renormalized interaction potential. Expectation values in the matrix model ()3.2fij) are 
defined in the usual way as 

0(M,Mt)\ /"DMDMt e-^^'^(*^^^W*^'*')) 0(M,Mt) , (3.28) 

/ E Zn{E) J 

such that they are well-defined and finite in the limit A^ — > oo. For later convenience, we 
have introduced in ()3.26|) a more general complex matrix model depending on a generic 
N X N, Hermitian external field E, which after calculation should be set equal to the 
diagonal matrix 

Sin. = I + /i' 1 (3.29) 



19 



appropriate to the original field theory, with 1 < i,m < N. This provides a concrete, non- 
perturbative definition of the noncommutative quantum field theory. Of course afterwards 
the limit in which the ultraviolet cutoff ()3.24|) is removed should also be taken. 

In what follows we will not be very precise about the convergence of the functional 
integrals (|3.26|) and (|3.28|) . In the regularized theory, the convergence of the traces is ob- 
vious. This is very similar to the situation in conventional quantum field theory, whereby 
one usually has an action functional which is only well-defined on fields which are suf- 
ficiently smooth or of Schwartz-type, while at the quantum level the functional integral 
also takes into account less well-behaved functions. At the algebraic level, the large 
limit should be rigorously defined by using the appropriate analog of that described for 
the noncommutative torus in |39j . 

While the above prescription is not the only possible definition of the quantum the- 
ory, we will see that it leads to a physically sensible and well-defined formulation of the 
quantum field theory, which is moreover exactly solvable. For this, it is crucial that an 
overall factor of has scaled out in front of the action in ()3.26p . With this convention, 
the correct scaling of the free energy which is finite in the limit A^ — > oo is given by 

^NiE) = In Zn{E) . (3.30) 

From a field theoretical point of view, this follows from the fact that the free energy should 
be proportional to the number of one-particle degrees of freedom. In a free scalar quantum 
field theory, this would be given by Tr ln(— 9^ -|- /i^) oc V, where V is the spatial area 
and A^ the area in momentum space. By using ()3.24|) . this quantity is proportional to A^^ 
in the large A^ limit. We will examine the issue of other possible scaling definitions of the 
matrix model representation in Section [7| 

It is important to note here that, unlike the conventional large 6 scalar field theories 
studied in ^26j where derivative terms in the action contribute only at leading order in 
j^, the model we solve here retains its kinetic energy terms exactly at A^ = oo. The 
exact solvability of the model is not ruined in this case because of the persistent GL(oo) 
symmetry, contrary to the cases of where the kinetic energy breaks this symmetry. In 
fact, already at the classical level one can see some important differences. The classical 
solutions of the 6 = oo models studied in [211 sue given by Hermitian matrices X obeying 
the matrix equation V'{X) = 0. They are of the form 

Ar = ^/i^0£,£, (3.31) 

e 

where fig are the distinct real critical points of the potential V{w). The corresponding 
field configurations correspond to the standard 9 = oo CMS solitons 130] and possess 
an infinite degeneracy under [/(oo) rotations of the matrices X, leading to a non-trivial 
soliton moduli space which is an infinite-dimensional Grassmannian jH]. On the other 
hand, in the present case the classical equations of motion are 

S + V'{X) = , (3.32) 

where X = M^M. The de generacy of the soliton solutions (and the existence of the 
ensuing moduli space) is now lifted. The solutions of (|3.32|) are still of the form (|3.31|) . 
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but now the real numbers /x^ satisfy the equation £u + V'{iii) = 0. For any polynomial 
potential V, this latter equation at large i generically has only a single solution, since 
— > cxD as £ ^ cxD and V'{w) ~ gdegv W^"^^^'^ for \w\ — >• oo. Because of the kinetic 
energy term there exist only finitely many real soliton solutions. As we will see in the 
following, this fact also makes a big difference at the quantum level, wherein the saddle- 
point solutions may be regarded as perturbations around the GMS soliton solutions [20] 5 
and the functional integral as an integral over the soliton moduli space. Indeed, the 
expansion of functions on with respect to the Landau basis may be thought of as a 
"soliton" expansion of fields on the noncommutative plane in terms of projection operators 
and partial isometrics of the Heisenberg algebra. 

As we will demonstrate in the next section, the crucial simplification which occurs at 
the special point i? = 6'"^ is that the partition function and observables of the matrix 
model depend only on the eigenvalues of the external field E. This enables an explicit 
solution of the matrix model. In Section [SI we will then see that the computations of 
the Green's functions fj2.12|) boil down to extracting the matrix model solution, and 
sums over Landau levels of certain combinations of the Landau wavefunctions which give 
the spacetime dependence of the correlation functions. This will then enable a non- 
perturbative analysis of the possible scaling regimes in which the cutoff ()3.24|) can be 
removed. 



4 Exact Solution 

We will now show how to solve the quantum field theory at the special point in parameter 
space where B6 = 1. At this point we can use the correspondence with the large N limit 
of the complex one- matrix model to obtain an exact, non-perturbative solution of the 
noncommutative field theory. We will begin by showing how to solve the general external 
field problem defined by (|3.26p at finite A^, thereby establishing the formal integrability 
of the model. Then we shall demonstrate how to obtain the solution of ()3.26p in the large 
A^ limit, and from this extract the solution of the noncommutative field theory defined 
by the action ^H^ . 

4.1 Integrability 

We will first derive a very simple determinant formula for the partition function ()3.26|) at 
finite A^. While this expression will not be particularly useful for extracting the solution 
of the corresponding continuum noncommutative field theory, it will formally prove that 
the model is exactly solvable and provide insights into the potential extensions to more 
general exactly solvable noncommutative quantum field theories. For this, we will consider 
the more general class of noncommutative field theories with actions (j2.H) . The regulated 
field theory is obtained through a modification of the matrix model (|3.26p defined by the 
partition function 

Zn{e,E)= I DMDMt e-^Tr(M£Mt+MtEM+y(MtM)) ^ 
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where the external fields E and E have respective eigenvalues Xg and Am, i,m = 1, . . . ,N. 
The additional external field E explicitly breaks the U{N) symmetry of the original model. 
Remarkably, this extended matrix model is nevertheless still exactly solvable. For this, 
we use the polar coordinate decomposition of a generic N x N complex matrix M which 
enables us to write 

M = f/tdiag(ri,...,rjv)F, (4.2) 

where U and V are N x N unitary matrices, and ri G [0, oo). This transformation 
introduces a Jacobian 

N 

DM DM^ = [dU] [dV] Y[dhi A^[/i]^ , (4.3) 

i=i 

where [dU] denotes the invariant Haar measure on the unitary group U{N), = rj, and 

l<£<m<Af 

is the Vandermonde determinant. Here and in the following we will not write irrelevant 
numerical constants (depending only on A^) explicitly. Upon substituting ()4.2|) and ()4.3p 
into ()4.H) . we observe that the integrations over U,V & U{N) decouple over the two types 
of external field terms in the integrand. With H = diag(/ii, . . . , hjsf), the integration over 
the angular degree of freedom U may be done by using the Harish-Chandra-Itzykson- 
Zuber formula |121 ES] 

[dU] ^-N^riEU^HU) ^ 1 l^-NX,K^^ (4 5) 

An analogous formula holds for the integration over V with E and Am- By expanding the 
determinants into sums over permutations of iV objects and using the antisymmetry of 
the Vandermonde determinants ()4.4|) . the partition function (14.11) may thereby be written 
as the eigenvalue model 

Z^[e,E]=Y, ^"^^ . J] / dh, . (4.6) 

^esjv ^iv[A] Ajv[A] J 



The eigenvalue integral ()4.6|) illustrates the necessity of the scaling ^ ~ in ()3.24|) in 
the large A^ limit. It is the only scaling that correctly takes into account the classical and 
quantum fiuctuation terms in ()4.6p . such that the effective action for the A^ eigenvalues 
hi is of order A^^. It may be expressed in the compact form 



E,E]=——^—-^ det ^ (/(A, + Am) ) , (4.7) 



Ajv[A]A7v[A] l<tm<N 

where we have introduced the special function 

oo 

'N{v{h)+Xh) 



/(A) = / dh e-^^^W+A^; _ (4,^ 
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Note that here V can be any function such that the integral in ()4.8|) makes sense. 

One way to understand the exact solvabihty of this model is to note that the partition 
function (j4.7|) provides a non-trivial solution of the two-dimensional Toda lattice hierar- 
chy [44].^ For any polynomial potential V, we may write the determinant formula (j4.7|) 
formally as 

Zj,(e,e] = det I <(> I ^~/fiz-^ + z-^) 

V ' / o<^,m<Jv-i \ J 2-niz J 27Tiz ^ ^ 

\ Z=0 2=0 

oo \ 

X J]e^(*^^'+*^^"M , (4.9) 

k=l / 

where we have defined the time variables 

N N 

e=i m=l 
with k > 1. This formula is derived by using the Cauchy identity to get 

oo N ^ 

M fui-zXf)(i-zXf) ' 



k=l 



and expanding the integration contours in ()4.9|) to catch the residues of the poles at 
z = 1/Xe and z = 1/Xm- This makes Zn{E,E) a tau-function r[t,t] of the integrable 
two-dimensional Toda lattice hierarchy. Its large N limit may thereby be determined 
from the equations of the dispersionless Toda hierarchy , although we shall not pursue 
this interesting line of approach to solving the noncommutative field theory ()2.1|) in this 
paper. A more transparent indication of the exact solvability of the matrix model (j4.H) 
is described in Appendix C. 

Let us now go back to the original model ()3.26|) . which is the E = limit of ()4.ip and 
hence should correspond to a particular reduction of the Toda lattice hierarchy above. 
In this limit, the expression ()4.7|) for the matrix integral is indeterminate. One can take 
care of this by regularizing the E limit in any way that removes the degeneracy, and 
applying rHopital's rule to (j4.7p . In this way we get 



where 



V 7 (_iv)"^ dX"" J ^ ' 



■^The analysis in the remainder of this subsection is geared at integrabihty features. It is not required 
in the remainder of this paper and may be skipped. The important result which we will need is the 
derivation of 14.25|) . An alternative derivation of this equation is given in Appendix C. 
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Thus once the function ()4.8|) is known exphcitly, then so is the partition function ()3.26|) 
for all N. The determinant formula ()4.12j) implies that the partition function Zn{E) 
is a tau-function r[t] for the integrable KP hierarchy |46j in the times tk of ()4.10j) with 
one-particle wavefunctions f^^\X~^), f^^\\~^), f^'^\X~^), . . .. This follows from the fact 
that, for any polynomial potential V, the functions ()4.13|) have the correct asymptotics 

^£-1 (^-1) ^ ^£-1 (^^ ^ (^-1) j as A ^ OO (4.14) 

such that (|4.12j) is a KP tau-function. Note while generic external field matrix models 
have other properties and correspond to particular reductions of the KP hierarchy the 
partition function ()3.26|) is exactly a KP tau-function. Its H^i+oo symmetry then manifests 
itself as the symmetry of the noncommutative quantum field theory under area-preserving 
diffeomorphisms which was described in Section 

We can also see this connection at the level of the equations of motion. The Schwinger- 
Dyson equations for the matrix model ()3.26p follow from the identities 

/ DM DMt y^^fu,^, ^-NT,{MEMf+ViMm))\ ^ Q ^ 
J ^ oMirn ^ ' 

which give rise to a system of differential equations for the partition function ()4.12j) . 



1 ^ 

N ^ 



d 



^ N dE ' 
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Z^iE) = . (4.16) 



Both the integration measure and action in ()3.26j) are invariant under arbitrary unitary 
transformations of E. The partition function Z]sf{E) = Z[X\ thereby depends only on the 
eigenvalues Xi of the external field E and is symmetric under permutation of them (This 
is manifest in the determinant form ()4.12|) ). Thus only N of the N"^ Schwinger-Dyson 
equations ()4.16p are independent. Writing the diagonalization explicitly as 

E =WUiag(Ai,...,A^) W (4.17) 

with U unitary, by using the chain rule the derivatives appearing in ()4.16|) may be rewritten 
via second order perturbation theory according to 

Let us now specialize to our main model of interest, the noncommutative complex $^ field 
theory defined by the potential ()2.9j) . In this case ()4.8j) is a transcendental error function 
that can be written in the form ^42j 




where a(C) is a monotonic function that increases from a(0) = ^ to a{oo) = 2. It is 
important to note though that the integrability property holds for any potential V. The 
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exact solvability of the noncommutative field theory is a consequence of its underlying 
matrix representation, rather than of a special choice of interaction. 

In this case, written in terms of the eigenvalues, the N independent equations of 
motion read 



_g__^ 9 ^ 1 ( d d_\ 1_ d 



Z[\] = . (4.20) 



Contracting (j4.2Up with A^'^' and collecting coefficients of the terms Xj'/N'^ leads after 
some algebra to a set of differential constraint equations in terms of the time variables tk 
of KW} of the form 



oo 



k=2 



k \ r\ r\ O O O 

O O sr^ O O O 



J2^k+k'-2 U[i]Z[A] = 0, (4.21) 

where 

•"A^W =9y2ir 7^ + SZ'^^" 7^ + T^r + ^ - 24,2 . (4.22) 

^ dtn dtk-n ^ dtn <9tfc dtk-2 

n=l n=l 

These equations coincide with the Virasoro constraints of the KP hierarchy ^H]- However, 
while all of this provides a nice characterization of the exact solvability of the underlying 
matrix model, it is not very convenient for dealing explicitly with its large N limit, which 
is what is needed for the precise connection with the original noncommutative field theory. 
The way to tackle the — > oo limit is the subject of the next subsection. 



4.2 The Master Equation 

To deal with the solution of ()4.20|) in the large limit, we define the function 

and introduce the spectral density of eigenvalues of the external matrix E, 

1 1 ^ 

p(A) = ^ Tr 5(A - E) = - ^ 6{X - A,) . (4.24) 

e=i 

At = oo, the eigenvalues Xi are assumed to be of order 1, and (j4.24p becomes a 
continuous function of A with support on some finite interval [ai,a2] on the real line. 
Setting Am = A in fl4.2H|l . a simple power counting argument shows that the derivative 
dW/dX which arises in the first term of ()4.20|) is of order jj: and so can be dropped at 
N = oo. Then, the differential equations ()4.2Up become 

- ~g W\X) -~g JdX' p(AO ^^^'^^ ~ ^^^^ + A W^(A) + 1 = 0, Xe[a„a,]. (4.25) 
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It is straightforward to obtain a perturbative solution to the non-hnear integral equa- 
tion ()4.25|) by writing a power series expansion for the function ()4.23p . 



A 

k-l 



1 

1^(0 (A) vr(^-'-i)(A) 



iy(A;^) = ^^'^W^«(A). (4.26) 

k=0 

By substituting this expansion into ()4.25|) . we can then compute the iterative solution 

1=0 

+ p(A') ^ ^-f ^ (4.27) 

ai 

with A G [ai, 02] and k > 1. The free solution W^^\X) of course follows directly from the 
Gaussian, g = partition function 

ZP{E) = e-^TrinE _ (4 28) 

An alternative interpretation of the function ()4.23|) as the resolvent of a related Hermitian 
one-matrix model is given in Appendix C. 

The exact, non-perturbative solution of ()4.25j) is worked out in Appendix D. For the 
logarithmic derivative of the partition function, we find 

yyi^X ) = _ V (^"i - bi){Xm - h) 



N 



1 a/ {Xj - bi){Xi - &2) - a/ (An, - bi){Xm - ^2) . . 

2iV ^ (A^ - A,) ^{Xe- b,){X,-b2) 

The parameters bi and 62 of this rational solution are determined by the non-linear con- 
straints 

N 

N t[^{X,-b,){X,-b^) 
V(A^-6i)(A^-62) 



4.3 The Exact Vacuum Amplitude 

The solution of the previous subsection represents the large limit of the general external 
field matrix model ()3.26|) . For the noncommutative field theory of interest, we should now 
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set the external field E equal to ()3.29p . For this, we shift the eigenvalues 
in (I4.2fjl4.3()l) . and set 



r,2 



X„ 



m 



Am A'' 



(4.31) 



where we ignore the zero-point energy shift which vanishes in the large N limit. We then 
need the function W{\) in the limit N ^ oo, m ^ oo with ^ e [0, 1] fixed. The spectral 
density 



P(A) 



167r 



(4.32) 



is flat, and it is supported on A G [0, IGtt]. The resulting integrals appearing in ()4.29l4.3Up 
can be done explicitly, and after some algebra we arrive finally at 



W^b(A) 



A + /i2 v^(A-6i)(A-62) 



+ 



IGtt 



In 



v/&i(A-&2)- V^2(A-6i) 



V(&i - 167r)(A - 62) - V(&2 - 167r)(A - 61) 



(4.33) 



where the parameters 61 and 62 are determined through 
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61 + 62 + 2/i^ = ^ In 
47r 



(61 + &2)' 



9_ 
27r 



— IGtt — \/h2 — 167r ' 

- V(&i- 167r)(62- 167r)) - 3/2^ - 8^ 



(4.34) 



The free energy — In Zb of the noncommutative quantum field theory, representing the 
exact (connected) vacuum amplitude, can now be obtained by integrating the function 
()4.33|) over A G [0, IGvr]. Using the permutation symmetry of the finite partition 
function Z[A] in the eigenvalues A^ and the boundary conditions ()4.34p . we arrive after 
some algebra at 



InZh 



hm — — 

1 

32^ 



N 



\/hh2 + ^{h - 167r)(62 - 167r) 
X ln4 [yW- ^/hi^ (^^Jhl-lQ^l- ^^2 - 167r) 
- ^/hh2 1116162 - a/(6i - 167r)(62 - IGtt) ln(6i - 167r)(62 - lOvr) 
+ ^ (Afi^ + 64 + 61 + 62-2 v/(6i - 167r)(62 



167r) 



4^ 



(61 + 62) (61 + 62)' - 3/iM + (6 



r.2 



2(61 + 62) +4/i 



(4.35) 



The loop amplitude (|4.33|) is a multi- valued analytic function on the complex A-plane with 
a square root branch cut along the interval [61, 62]. From the constraint equations (j4.34j) . 
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it is straightforward to see that the two branch points 61 and 62 are always complex- 
valued. Moreover, the one-cut solution ()4.33p is an analytic function of the parameters ft'^ 
and g. Thus the matrix model exhibits no non-trivial critical behaviour as the coupling 
constant g is varied, and there is no scaling limit at large with an approach to a phase 
transition. This is completely consistent with the fact that the regulated noncommutative 
field theory already describes a continuum model. The problem of removing the ultraviolet 
regularization (j3.24|) will be analysed in the next section. 



5 Correlation Functions 

In this section we shall investigate the spacetime dependence of the Green's functions 
(j2.12|) using the exact, non-perturbative solution of the matrix model that we obtained 
in the previous section. As an illustration of the power of working in the Landau basis 
for the expansion of noncommutative fields, and also as a warmup to some of the general 
effects that we will see, we shall first investigate a novel property of these correlators 
in a particular truncation of the field theory which has an extended GL {00) x GL(oo) 
symmetry. Then we will proceed to obtain exact expressions, valid to all loop orders, for 
the connected Green's functions of the noncommutative quantum field theory. This will 
enable us to describe explicitly the non-perturbative scaling limits of the model. 



5.1 Green's Functions in Static and Free Limits 

An important feature in the analysis of the Green's functions fl2.12|) is that we use ordinary 
products of the fields. A natural question which then arises is as to what the difference 
is between these functions and those defined using star-products, as may seem more 
appropriate to a noncommutative field theory with a star-unitary symmetry. To get a 
fiavour for the answer to this question, we will find an L^-integration kernel . . . ,xl) 

with the property that 

d'^x fi-k- ■ - -k fL{x) = j d^xi ■ j d^XL ^L(a;i, . . . /i(a;i) ■ ■ ■/^(xl) (5.1) 

for any collection of L fields /i, . . . , f^. For L = 2 this kernel was computed in ^2]- The 
identity (j5.1|) thereby produces a convolution formula 

Yl d^xi d^yi g2r{xi,yi,...,Xr,yr) G^^\xi,...,Xr]yi,...,yr) 
1=1 

= j d^x {{<^U<^)[{x)^ (5.2) 

for the integrated Green's functions of the noncommutative quantum field theory in terms 
of correlators of composite operators. 

The kernel Ql bears a remarkable relation to the matrix model constructed in this 
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paper. Defining 5y{x) = 5^'^\x — y), from tlie definition ()5.1|) it follows that 

^^(xi, . . . = y d^a; 5^., ^ ■ ■ ■ * 5^.^ (x) . (5.3) 
Using completeness of the Landau wavefunctions on to write 

oo 

e,m=i 

and the projector relations ()3.12|1 . we may compute the star-products appearing in ()5.3j) 
according to 

^ oo oo 

Sxi ^S^ii^) = nr-?\ 4>i,rrXA V'</'<?,fc(a;2) (^fc,m(a:i) • (5.5) 
v47rt^ f- — \ ~, 

By doing this an additional L — 2 times in ()5.3j) . and then integrating over x using the 
orthonormality relation ()3.1()|1 . we arrive at a trace formula for products of the Landau 
wavefunctions, 

L-2 



\ n—z oo 

/— - j 0mi,m2(2^L) 0m2,m3(a;L-l) ■ ■ ■ 0m.i,,mi(a^l) 

' mi,. ..,1711^=1 

TT^(^(l){xL)(t){xL-i)- ■■(t>{xi)^ . (5.6) 



mi,...,m£=l 

The physical significance of this formula is clear. Before regularization and at 5 = 6^^, the 
— ^ oo limit of the matrix model action (j3.18|) is an invariant function of the combination 
M'^M and thereby possesses an extended GL (oo) x GL{oo) symmetry M y-^ U ■ M ■ V~^, 
M"^ I— V ■ M'^ ■ U~^. In this limit, the kinetic energy completely drops out and the 
field theory possesses only static configurations. Because of this huge symmetry, the only 
non-vanishing correlators of the large N complex matrix model are those which depend 
solely on traces of powers {M'^ My . Using the field expansions ()3.13p . it follows then that 
the right-hand side of (|5.6p gives the exact spacetime dependence of an L-point Green's 
function (|2.12j) in this limit, with L = 2r. In other words, the solution to the star-product 
relation (j5.1|) is determined by the bare Green's functions in the limits of interest here, 
obtained essentially by dropping the kinetic energy terms in (|2.H) . i.e. a = a = 0. 

The integration kernel ()5.6|) is computed explicitly for all L in Appendix E. For the 
cases L = 2r which are relevant for the complex scalar field theory, we have 

Q2r{Xi, . . . , X2r 




X exp 



- Y^{-iy-^{J-I-r)xj-Bxj] . (5.7) 



r 

KJ 



Note that this result holds for any interaction potential V. The full Green's functions 
(|2.12p in the given limit are obtained by taking the sum of (|5.7p over all permutations 
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of the coordinates xi, . . . ,X2r,^ and then multiplying it by matrix model correlators of 
Tr (M'l'M )*" which carry the detailed information about the potential ()2.4|) of the quantum 
field theory. These matrix averages are well-known to be calculable through a generating 
function which is defined by the contour integration [15] 

^E;^(^(^^'^0'L„. (5.8) 

r=0 

where the normalized expectation values on the right-hand side of (j5.8p are defined by 
()3.28|) in the large limit, and the parameter a is determined in terms of the rescaled 
coupling constants of the potential ()3.27p as the real positive solution of the equation^ 



1 

2 



dw wV'{w) fz 



[—a,a\ 



2 , 

k>2 k 



[{k-l)\_ 



4fc 



The connected part of the 2r-point Green's function is of order 1/N^~'^ in the large 
limit, due to the usual factorization of correlators in the matrix model at = cx3. We will 
elucidate this point further in Section 15.31 Notice also that the large 6 bare propagator 

g2{x,y) = 5^^\x-y) (5.10) 

is ultra-local due to the persistent GL{oo) x GL{oo) symmetry.® This is no longer true, 
however, for the higher order Green's functions, which respect the magnetic translation 
symmetry described in Section] 



Now let us examine the opposite extreme where the interaction potential is turned off, 
and compare these exact results with the free propagator 

C,ix,y) = {x\:^-^\y) (5.11) 



"'Note that (|5.7|) is not real-valued, because even for real functions / and g, f-kg — g-kf^^f-kgm 
general. It becomes real after symmetrization over the spacetime coordinates. 

^This result is derived by introducing the resolvent function R{z) = -j^ (Tr ■^2rrfjTM) e=o' "^^lose large z 
expansion coincides with the right-hand side of H5.8|l . One can then write an appropriate Schwinger-Dyson 
equation analogous to (|4.15|) for this correlator, which at = oo takes the form 0S] 



/dw wV'(w) , , „ / 1 

— „ ^ „^ R{w) = R^{z) + 0[—r 
47ri z2 _ ^2 V y \ J yj^2 



where the contour C encloses the singularities of R{z), but not the point w — z, with counterclockwise 
orientation in the complex z-plane. Assuming that the only singularity of R{z) is a branch cut across a 
single connected interval [—a, a] C M, the solution of this equation coincides with the left-hand side of 
(|5.8(l . Then (|5.9(l follows by imposing the asymptotic boundary condition R{z) — + 0{z~'^) as z ^ oo 
on this solution. 

^The form (|5.1U|) follows immediately from H5.6|l by completeness of the Landau wavefunctions, 
or alternatively from the definition (|5.1|l by the well-known star- product identity J d^x f * f'{x) = 
Jd'x fix) fix). 



30 



of a charged scalar particle of mass /i in the constant electromagnetic background. The 
computation of ()5.11|) is presented in Appendix F and one finds 



1 1^ ..12 /na ;„ D.. / - e / U \ / 



Air 







.^u^ + \x -y\'^u/e \u+\x-y\ye 



(5.12) 



In particular, by setting /i^ = in ()5.12p . the integral can be done explicitly and yields 
the free massless propagator 

Coix,y) = ^ e-'^-^'y Ko{\x-y\y2e) , (5.13) 

where Kq is the modified Bessel function of the second kind of order 0. From ()5.12|) we 
see that the ultraviolet problem is unchanged by the presence of the electromagnetic field 
B = 1/6, i.e. the singularity at x = y is the same as that for B = 0, while from (j5.13|) we 
see that the magnetic field acts like a mass and cures the infrared problem in the massless 
case. In particular, as shown in Appendix F, the 6 ^ oo limit of the free propagator yields 
an ultra-local spacetime dependence only for massive fields /i^ > 0, while the result ()5.1()j] 
is valid for all values of the coupling constants of the field theory. In the massless case 
()5.13|) . the standard logarithmic ultraviolet behaviour appears only at very short length 
scales |a; — ?/| <C ^/6. As we will see more of in the following, these properties are the 
remnants of the UV/IR mixing phenomenon within the present class of noncommutative 
quantum field theories. 



5.2 The Exact Propagator 

In this subsection we will compute the exact two-point function of the quantum field 
theory, which by using the field expansions ()3.13|) can be written formally as 

oo 

G{x,y) = (^<l>\x)^y))=A7fd ^ (m^,, M,,^,)^^^ 0,,^(x) ,,(y) . (5.14) 

e,m,e',m'=i 

To make sense of ()5.14j) . we restrict the sums to < N (the sums over m and m' 
need no regularization) and compute the matrix correlator in the 't Hooft limit of the 
matrix model. The action and measure in the functional integral ()3.28p are invariant 
under unitary transformations ()3.21|) with W = U^^. We can thereby explicitly make 
this transformation in the matrix integral and then integrate over the unitary group. At 
finite the integral we need is jlH] 

/ W]UlU,,, = ^5uy5u' , (5.15) 

{/(TV) 

which can be derived by using the invariance properties of the Haar measure. As a 
consequence, the matrix expectation values appearing in ()5.14|1 are given by 



(5.16) 



E=S 



31 



Using ()4.23p and ()4.3H) . the Green's function ()5.14|) is therefore given by the large N hmit 
of 



GN{x,y) = -- J2w{\i)leix,y) , (5.17) 
where we have introduced the sum over Landau levels 

oo 

Ji{x,y) = 4:719 ^ 4>mA^) (Pi,m{y) ■ (5.18) 

m=l 



The function (j5.18p is computed in Appendix G. It follows from ()G.4|) that, at finite N, 
the propagator ()5.17p respects the magnetic translation invariance described in Section lT^ 
As discussed there, however, this will no longer be true in the desired limit — > oo with 
6/N finite. The two-point function ()5.17|) is to be evaluated in the limit A^ — oo, £ ^ oo 
with G [0, 1] fixed. The function Vr(A^) in this scaling limit is given by ()4.33p . while 
the sum over Landau wavefunctions ()5.18|1 is computed in Appendix G with the result 

7^(x,t/) = 4Jo(^Av^|x-i/|) , (5.19) 

where Jo is the Bessel function of the first kind of order 0. In the large A^ limit, the 
function ()5.17|) thereby becomes 

Wtt 



G{x,y) = - J^W^iX)Jo(^Ay/X\x-y\) 





40FA 



/ ^^W^ip'/^') Jo{p\x-y\) , (5.20) 







which by using the angular integral representation of the Bessel function can be written 
in the form 

G{x, y) = -Y2 J ^2^^ (pVA^) e . (5.21) 

|p|<40FA 

The result (15.211) illustrates, along with the Green's functions of the previous subsection, 
a remarkable property of the large 6 limit of the quantum field theory described in Sec- 
tion EUl In this limit, the underlying spacetime is expected to disappear or to degenerate 
and all spacetime symmetries to be maximally violated. Instead, the Green's functions 
are both rotationally and translationally invariant. This is characteristic of a "fuzzy" 
regularization scheme, which typically preserves all spacetime symmetries of the original 
continuum field theory (in contrast to lattice regularization). The remnants of UV/IR 
mixing here appear in the far infrared at |a; — 2/| ~ y/O. In deriving ()5.21|1 . these distances 
have been effectively scaled out, so that all results here are valid at — y| ^ VO- Thus the 
spacetime picture which emerges in the limit described in Section EiSl is not pathological. 
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From ()5.2H) it also follows that A is clearly the ultraviolet cutoff in the quantum field 
theory. More precisely, the quantity 4-(/7rA is a sharp cutoff in momentum space. Thus 
the large 9 limit defines a rotationally and translationally invariant field theory with a 
finite cutoff. Moreover, the exact solution of the matrix model Schwinger-Dyson equations 
has a direct physical meaning as the exact propagator of the quantum field theory. From 
(j5.21|) we see that the function (j4.33|) is essentially the Fourier transform of the exact 
two-point function, 

G(p) = -^l^b(pVA') (5.22) 

with IpI < 4 -^/tF A. The crucial issue now is whether or not the momentum cutoff can be 
removed, i.e. if the limit A — > cxd can be taken in the interacting quantum field theory. It 
is straightforward to see that the two-point function is not finite in the naive scaling limit 
fl = fi/A, g = g/ h?. This can be seen directly in perturbation theory. Using the iterative 
solution (j4.27p and (j5.22j) . we can easily determine the propagator up to two- loop order 
in the limit A — oo as 



G{p) 



^ ln(167rAV/i2) g 



g^ /ln(f67rAV/i2) A^ 3^^ i^^^g^ ^2/^2^ 



+ 



167r/i2(p2 +^2)2 \^ 2/iV 167r(p2 + ^2)4 

g"" /ln2(I67rAV/i') , , 2g^ 

f Gtt {jP + n 



2^3 



V 16,1 ^ ^^)^ (rf + i^^f^ ' ■ ' ' 



The first term in ()5.23|) is the free propagator which as expected is finite as the momentum 
cutoff is removed. The second term recovers the usual, one-loop logarithmic ultraviolet 
divergence of theory in two dimensions which is generated by the planar bubble dia- 
gram. However, at one-loop order there is an additional quadratic ultraviolet divergence, 
which may be attributed to the non-planar contribution in this case. The divergences are 
even worse at higher loop orders. For instance, the usual two-loop ultraviolet divergence 
generated by the planar sunset diagram is accompanied by terms in (j5.23p whose degree 
of divergence is different than that in the usual scalar field theory. These additional di- 
vergences in A are nothing but divergences in the summations over degenerate Landau 
levels in the matrix model, whose degree grows with the order of perturbation theory. 

Another way to try to get rid of the cutoff is to adjust the coupling constants non- 
trivially. For this, one needs to find a critical regime in which there is another scale that 
is smaller than A and the couplings can be thereby adjusted to achieve a separation of 
scales, such that non-trivial physics remains for momentum modes with \p\ <^ A. With 
the exact solution at hand, it is not difficult to show that this is in fact not possible. 
For this, we note that from ()5.22|) it follows that the behaviour of the propagator at 
large distances (much larger than the inverse momentum cutoff A~^) is determined by the 
singularities of the loop amplitude W{X), i.e. its two complex branch points A = &i,62- 
As a consequence, we can write down the asymptotic behaviour 

^(a;,?/) ~ e-l^-^l/« , \x ~ y\ ^ A-\ (5.24) 

where the correlation length ^ = 1/Impo is determined by the condition that the complex 
number z = Pg/A^ coincides with one of the two branch points bi or 62- From (j4.34j) one 
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can show 

(Im zf -2iRez)^ = blbl-^ (^61 + 62) % 2^ . (5.25) 

From this inequahty one may then show that for any choice of parameters, the correlation 
length ^ is bounded as 




(5.26) 



From ()5.2(ij) it follows that the only way to keep the correlation length finite while 
sending the cutoff to infinity is to send the coupling constant g to zero. Careful inspection 
of the function ()4.HH|1 shows that the only scaling limit which makes sense is ^ ~ 
which is also suggested by the inequahty (j5.2(ij] . Thus we define 



4 

m 



9 = (5-27) 

with m a dynamically generated mass scale. From the constraint equations (j4.34j) we 
further find that 6162 = 2^ and also fei + 62 = 0{g), which drops out in the scaling limit 
described above. The renormalized two-point function is thereby given from (|4.33p and 
(KT]^ as 



SO that the additional limit A ^ 00 makes the spectral density of the external field invisible 
in the loop equation ()4.25|) . The meaning of the scaling limit in ()5.28p is easy to under- 
stand from perturbation theory. Comparing the large momentum expansion of ()5.28p 
with ()5.23|) . it is evident that taking this scaling limit amounts to a non-perturbative re- 
summation of the leading power divergences in perturbation theory. As these divergences 
come from the Landau level degeneracies and are thereby an artifact of the particular 
regularization used, such a resummation procedure does not lead to a meaningful inter- 
acting quantum field theory. This will be confirmed in the next subsection, wherein we 
show that all higher-point Green's functions are trivial. 

We conclude that the present model is not a renormalizable, interacting quantum 
field theory, because its coupling constant flows only to the trivial Gaussian fixed point. 
Within our present formulation of the model, there are a priori two scales at hand. One 
is the size of the orbit in the lowest Landau level, which is the infrared cutoff and 

which we have taken to be infinitely large. The other is the size ^JN/ B of the largest 
Landau level that we allow for, which is the ultraviolet cutoff. What we have found above 
is that there is no intermediate scale in between the two, no matter what one does with 
the couplings of the model. The fields are thereby correlated on the scale of the cutoff, 
and the field theory is not renormalizable in the standard sense. Of course, one needs to 
analyse the possibility that other scalings, besides the standard 't Hooft limit which was 
described in Section may lead to non-trivial fixed points of the renormalization group 
flows. In Section [7| we shall address this important issue and describe some ways that a 
non-trivial interacting quantum field theory may be attained. 
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Nevertheless, the present model provides an interesting example of a noncommutative 
field theory which has a finite cutoff, and which is exactly solvable. In Section [7| we will 
offer an alternative physical interpretation of the regularization of Section lT^ in which A is 
simply the characteristic mass scale of the theory, and for which the correlation functions 
at finite A can be interpreted as those of the full quantum field theoretic limit. The Green's 
functions display some novel features which are not encountered in ordinary quantum 
field theory. For example, because the branch points bi and 62 are complex-valued, the 
propagator exhibits an unusual oscillatory behaviour on top of its exponential decay. 
This is similar to the Aharonov-Bohm phases, for charged particles in the background 
magnetic field B = 6~^, that have been observed on top of the usual area law behaviour 
of Wilson loops in two-dimensional noncommutative Yang-Mills theory [SD] . The origin 
of this behaviour can be understood from the bare Green's functions ()5.7p . and the free 
propagators ()5.12j) and ()5.13|) for bosons in a finite external magnetic field, where the 
Aharonov-Bohm phase factors x ■ By are explicit. Similar behaviours for the Green's 
functions of scalar field theory on more complicated noncommutative spacetimes have 
also been observed in 



5.3 Connected Correlators 

To substantiate the claim made in the previous subsection about the triviality of the 
interacting quantum field theory, in this subsection we show that the four-point Green's 
function is trivial, i.e. it factorizes into products of propagators in the scaling limit 
described above, and further argue that this is a generic feature of all higher order corre- 
lators ()2.12|) . For this, we expand the four-point function analogously to the propagator 
to formally get 



00 

£1 ,mi ,.../4,m,4=l 

X (j)ei,mi{x) 0^2/2(1/) 4>i3,ms{z) (prrn/iiw) ■ (5.29) 



As in the previous subsection, we make sense of ()5.29|) by restricting the sums over ij 
to 1, . . . ,N and compute the matrix model correlators in the 't Hooft limit. Again, we 
replace M ^ U ■ M, M'^ M'^ ■ W and integrate over all A^ x unitary matrices U. A 
straightforward computation using the unitary matrix integral formula 



[dU] Ul^i,^ U^^f^^ ^tih ^Uki — jy2 



Slie2 ^kik2 ^i-iiA ^k^ki + ^£1^4 '^fci/c4 ^£2^3 ^k2k3 



U{N) 



(^^iiei ^kik2 ^hia ^kzkA + ^tit2 ^kxkA ^Is^a ^fefcs j 

(5.30) 



yields 
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1 ^ 

£l,.../4 = l 



(5.31) 



where the notation (< — indicates to include the same terms but with the spacetime 
dependence {x, y); {z, w) replaced by (x, w); {z, y), and we have introduced the generalized 
Landau sums 

oo 

Ju^ix, y) = AttO E (j)e',m{y) ■ (5.32) 

m=l 

The diagonal elements of the functions ()5.32p of course coincide with ()5.18p . 7^(x, y) = 
7«(a;,2/). 

We can now compute the matrix model correlators appearing in ()5.3ip in terms of the 
free energy ()3.30|) as 

dEi^e-^ dEe^e^ OEi^i-^ dEe_^i^ 

Since the function En{E) depends only on the eigenvalues Xi of the external field E, we can 
compute its derivatives by replacing E with E + v/N and using Schrodinger perturbation 
theory to expand it about v = (see ()4.18|) ). One finds 

d^NjE) 1 

-Olm Oi'm' yV2[M, M' ) 77 Olm' Oe'm I -L — Oem 



(5.34) 

where W{Xi) is the loop function ()4.23|) and we have defined 

The matrix correlator ()5.33|) is thereby given as 

(MtM)^^^^ (MtM)^^^^)^ = 6,,,,6,,,,(wiX,,)WiXe,)-W,iX,,,Xis) 

1 . . A ^ W^(A.i)-W^(A£3) 

77 Of,!', Of^f„ \ i — Of 



N ^^^^ V ^^^V A^, - A,3 

(5.36) 
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The crucial point now is the behaviour of the sums over Landau eigenfunctions in 
fl5.32|) in the large N limit. In the previous subsection we saw that for i = i' this sum 
is of order 1. For £ ^ f, hy using the analysis of Appendix G it is possible to show that 
the sum is exponentially small in \i — i'\ for i — i' c:^ N. Putting everything together, the 
Green's function fl5.31|) at = oo becomes 

1 ^ 

i,Tn=l 

= G{x,y)G{z,w) + Gix,w)G{z,y) + o(^-^^ (5.37) 

with G{x,y) the two-point Green's function ()5.17|1 . The right-hand side of ()5.37j] is just 
the "trivial" connected part of the correlator up to order j^, and this proves that the 
connected four-point Green's function of the field theory is trivial in the large limit. 
From the analysis above we believe that this triviality is unavoidable in the 't Hooft limit, 
as it is simply a consequence of the well-known large A^ factorization in the matrix model, 
which implies that all connected matrix correlation functions vanish. In particular, it 
should be valid for arbitrary interaction potential V. We believe that this is a generic 
feature that will generalize to all higher-point functions of the model, and so the quantum 
field theory in the scaling limit obtained in the previous subsection is Gaussian, as all 
connected 2r-point Green's functions are of order l/N^"^. 

The computation above shows that the connected part of the four-point Green's func- 
tion, 

<^conn(a;, y; z, w) = G'^'^\x, y, z, w) - G{x, y) G{z, w) - w) G{z, y) , (5.38) 

is of order in the 't Hooft limit. However, it also shows that one can compute the 

quantity N Gtonn{x, y; z, w) explicitly in the large A^ limit, and everything needed to write 
down a closed formula for it is given above. The non- vanishing leading term of this and 
higher Green's functions may be of interest when applying this model as a generalized 
mean field theory of a local quantum field theory, as we discuss in Appendix A. It provides 
the analog of the random phase approximation in conventional mean field theory. A more 
detailed investigation of these higher-order Green's functions would be interesting but is 
left for future work. 



6 Higher Dimensional Generalizations 

A remarkable feature of the present class of field theories, and the formalism that we have 
developed to analyse them, is that everything we have said thus far generalizes to arbi- 
trary even dimensionality. Namely, the quantum field theory with action (j2.1Uj) is exactly 
solvable for any n. This is quite unlike the situation that one would expect in ordinary 
quantum field theory. In this section we will describe the higher-dimensional generaliza- 
tion of the exact solution that we have obtained. As we will see, the 2n-dimensional case 
can be mapped exactly onto the two-dimensional model (This is also observed in [T^ l^). 
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Since much of the machinery is the same as in the two-dimensional case, we will only 
highlight the essential changes which occur. 

Let us first fix some notation. We can rotate to a local coordinate system on R^" 
in which the non-degenerate antisymmetric matrix 9 = (^'-') assumes its canonical skew- 
diagonal form 



/ 







V 



\ 





/ 



(6.1) 



with non-zero skew-eigenvalues 6\ i = 1, . . . ,n. Similarly, at the special point B = (Bij) = 
6~^, which we assume throughout, the antisymmetric matrix B is in its canonical form 
with skew-eigenvalues Bi = (6'*)^^ Corresponding to the i^^ skew block, we introduce the 
complex coordinates 



(6.2) 



We will also use a similar notation for momentum space variables, Ki = {p2i-i,P2i) 



,n. 



6.1 Regularization 

The Hamiltonian — in the action 1)2.111 is the sum of n two-dimensional Landau Hamil- 
tonians in the magnetic field variables B^. Its eigenfunctions are therefore given as a 
product of two-dimensional Landau wavef unctions. 



l[^i,,m, (z\z^;B,) , (6.3) 



1=1 



where mi(2;*, z'^ ; Bi) are the two-dimensional Landau eigenfunctions of Sectional in the 
noncommutativity parameter 6'* = {Bi)~^, and the quantum numbers are 

£ = {ii, . . . ,in) , m = (mi, . . . ,m„) , fj, = 1, 2, . . . . (6.4) 

The functions ()6.3|) form a complete ortho normal basis in L^(M^"). From the correspond- 
ing two-dimensional results of Section El we obtain immediately the eigenvalue equation 



i=l 



- D^^,^ = Aj2Bi(ii- l) ^i,m , (6.5) 
and the star-product projector relation 

det(47r^)V4 



= :nr7JZ7rrui ^■m.,£' (6.6) 



38 



with Si^rn = rii^^i.mr All analogous expression holds for the Hamiltonian — D^. 
We can now diagonahze the action (|2.H) by the expansions 



(6.7) 



£,Tri6N' 



where s is a scahng parameter which we will determine below so as to get the appropriate 
scaling required for exact solvability of the model. We regard M = (Mg^) as a compact 
operator on the Hilbert space H = iS(N"), with trace Tr^(M) = Ylim^'m,m- We will use 

an obvious matrix notation with {M'^M)f^gi = ^rn,,^' ^^'-^ (^^)£,Tn, = 

Then the action 1)2.11) can be written in the matrix form 



where the external field is given by 



4 = 1 



2 1+" 



and we have defined the renormalized interaction potential 

Vo{w) = 



k>2 



det(47r^)('=-i)/2 k 



9k k 
w 



(6.8) 



(6.9) 



(6.10) 



There are many ways to now regularize the noncommutative field theory in the form 
(16.8)1 by mapping it onto a finite dimensional matrix model. As we will see later on, 
these regularizations are all equivalent in the limit where the size of the matrices becomes 
infinite. The simplest one, which we will call the "naive" regularization, is to restrict each 
set of Landau quantum numbers to a common finite range L < oo, 

ii,mi = 1, . . . ,L , i = l,...,n. (6.11) 

We can then map the integer vectors ()6.4)) bijectively to single integers i, m as 

n n 

i = Yl ^''^ (^4 - 1) + 1 , ^ = XI (mi - 1) + 1 . (6.12) 

1=1 i=l 

The precise form of the mapping 1)6.12)1 is not important here, as any other bijection 
on — >■ N would work equally well, for instance any lexicographic ordering t ^ I = 
lex(£i, ...,£„) and m \^ m = lex(mi, . . . , m„). With this mapping, we can then identify 
matrix elements as M^^n = M^^, and the regularization above amounts to restricting 
i,m = 1, . . . ,N with 

N = L'^ . (6.13) 

It should be stressed though that this is not the natural regularization. The one leading 
to a nice N oo limit is the isotropic regularization whereby the one-particle energies 
are bounded by a momentum scale A' as 

<A'^ (6.14) 
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The isotropic limit is essentially obtained by letting the skew-eigenvalues Bi of the mag- 
netic field all approach the same value. However, to get the correct scaling it is more 
convenient to proceed first via the naive regularization above. 

The regularized partition function is then given by the N x N matrix integral 

Zj,i£) = j DM DMt ^-^r{aMSMUaM^£M+V,{M^M)) ^ ^g^^g) 

similarly to that of Section f3.2[ To obtain a good large N limit, we now need an over- 
all factor of coming out in front of the trace in (j6.15p . This requires a relation 
det(i?)^/^"' L = 47riV = AttL"", which fixes the scaling parameter s introduced in (j6.7j) 

as 

s = L("-i)/2 det(47r^)i/^ . (6.16) 

The partition function ()6.15|) then assumes the form in ()3.26|) . with the renormalized 
interaction potential (j6.1(Jj) taking the form (j3.27|) . As expected from dimensional analysis, 
the rescaled coupling constants in (|3.27|) are now given by 



i2 

' ^'^ ^ ^2n-2(n-l)fc 

with 



The quantity ()6.18p is the natural ultraviolet cutoff in any dimension, as then the one- 
particle energies are bounded as 4 Biii < 4bLn = 167mA^ (in the isotropic case where 
all skew-eigenvalues Bi = h are the same). The one-particle energies are now given by 
^Lm = ^tmi with 



= 47r 

and 



i=i ^ ^ 



Bi 



det(5)i/ 



2n 



(6.19) 



(6.20) 



After shifting the eigenvalues ()6.19|) by the mass as in Section 14.31 they are scaled 
correctly so that the spectral density 



A^ 

is given by a finite Riemann sum in the large A^ limit, 

167r 



''W^^ E S(X-Xt) (6.21) 



"W^n /^^"^-E'"^.) ■ (6-22) 

i=l i \ i=l / 
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From ()6.22j) it follows that the isotropic regularization in ()6.14j) . i.e. ^^d-i < is simpler 
and more natural. In that case, rather than the limits Aj < IGvr in the integrals defining 
the density of states ()6.22j) . we have ^j&iAj < IGvr. We can then transform to spherical 
coordinates to explicitly compute the integral in ()6.22|) and obtain the large N eigenvalue 
distribution 

^n— 1 

p(A) = — ^ — ^ , (6.23) 

'^^ ' (16v^)" r(f) ^ ' 

with compact support A G [0, IGvr]. Here F is the Euler function which arises from 
integrating over the solid angle in n dimensions. Thus one of the main characterizations 
of the higher dimensional generalizations is that the spectral density is no longer fiat. 



6.2 The Two-Point Function 

Let us now consider the two-point function G'(x, y) in the case of the noncommutative field 
theory ()2.1()|1 . With the identifications £ = i, m = m, the matrix model representation 
f)6.7|6.8p is formally the same as in the two-dimensional case, and we thus find 

GN{x,y) = -^ J2 WiX,)^,ix,y) (6.24) 

with W{X£) the logarithmic derivative of the partition function (j6.15p in the 't Hooft 
scaling limit and 

7,(x,y)=det(47r^^)V2«L"-i J2 <Pi,U^) (pmAv) ■ (6-25) 

Since the Landau eigenfunctions in (j6.25|) are just the products (j6.3|) of two-dimensional 
Landau wavefunctions, the sum in (j6.25|) is a product of the n sums computed in Ap- 
pendix G and in the scaling limit whereby i/L is finite as L oo it gives the result 



n oo 



7,(a;, y) = det(47r0)i/2" L^-' J] ^ (z^ , ; 5,) (w' , ; B?j 

i=l m,i=l 

n 

= rA'''~'l[Jo[^/xJ^A\z"-w'\^ , (6.26) 



with = y^* ^ + iy'^\ = 2/^* ^ — iy'^\ and Aj = IGirii/L. In the limit L ^ oo, the 
Green's function fl6.24|) thereby becomes 



n 



G(x,i/) = -A2"-2 J] / -1 Jo^^/A^fo-A H^ K^M,, • (6.27) 

We now change variables AjfejA^ = Kj in ()6.27|) . use the fact that Hi^j — 1' ^^'^ then 
apply the angular integral representation of the Bessel function Jq to get 

G(-.!') = 4n / ^ /£e"^*'-'-w(Ef|) ■ (6-28) 
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We can thereby identify p2i-i = Ki cos Ti, p2i = Ki sin Ti as the two-dimensional mo- 
menta associated with the skew-block coordinates x^*), and write 

1 r d^"n 

y) = 'Y^ J J^n^ (pVA^) e -■(-^) , (6.29) 

{pi._i+pi,<16^A26,} 

which is the 2n-dimensional generalization of the formula ()5.2ip for the propagator. Thus 
the scaling limit used above always looks rotationally symmetric except at very small 
distances. However, this anisotropy is due to the anisotropic regularization that we have 
used, which is invariant only under an abelian SO (2)^ subgroup of S0{2n). Had we used 
the isotropic scaling described earlier we would have obtained a rotationally symmetric 
cutoff < IGttA^ in momentum space. No matter how one chooses the skew-eigenvalues 
Bi, in the large limit the anisotropy is totally washed out and one gets the same result. 
Thus, rather remarkably, we are left with an S'0(2n)-invariant Green's function, even 
though rotational symmetry is broken by the external fields. While other regularizations 
are possible as well, at finite cutoff A they only make a difference at small distances 
|x - yl < A"^ 

The Fourier transform W{X) of the two-point function satisfies the same loop equation 
fl4.25|l as in the two-dimensional case (before the mass shift), whose solution is given in 
fl4.29l4.30p . The only change is that the density of states p(A) is no longer flat, and 
with the isotropic regularization it is given by the simple explicit expression ()6.23|) . In 
particular, from the iterative solution ()4.27p it follows that the perturbative expansion of 
the propagator now contains the anticipated power law divergences of scalar field theory 
in 2n dimensions. Moreover, the analytic structure of the loop function iy(A) is fairly 
insensitive to the precise form of the spectral density, as long as its support lies in [0, oo). 
In particular, the singularities of W{X) are again cuts which sit in the complex plane, 
and it is possible to demonstrate once again that the field theory in any dimension is not 
renormalizable. In the limit A oo, the rotational symmetry breaking S0{2n) S0{2)^ 
is no longer washed away. 



7 Scaling Limits 

In the previous sections we have found that, in the scaling limit whereby the quantity 
fl3.24|) plays the role of an ultraviolet cutoff, the renormalized quantum field theory is 
trivial in the sense explained before. The purpose of this section is to point out that there 
exist other scaling limits of the noncommutative field theory which, while mathematically 
sharing the same properties as that studied already, have drastically different physical in- 
terpretations and have the potential of providing non-trivial connected Green's functions. 
At the same time, we will give nice physical characterizations of the large limit of the 
matrix model from a field theoretic point of view. We shall also describe extensions of 
our analysis that remove the Landau level degeneracies which spoilt the renormalizability 
of the interacting quantum field theory. 
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7.1 Scale Transformations 



The crux of the existence of different scahng hmits of the noncommutative quantum field 
theory resides in the special scaling property of the Landau wavefunctions, 

B) = A" (l)e,m (A X ; A"^ B) , (7.1) 

which easily follows from the generating function ()3.8|) . The physical interpretation of 
fl7.1|) is that a typical skew-eigenvalue b of B defines the magnetic length Zmag = ^/Vb. If 
the length scale is changed as Zmag A /mag, then the physics is unchanged provided we 
rescale the spacetime coordinates also as x — ^ Ax. The factor of A" in (|7.1|) then ensures 
the correct normalization of the Landau eigenf unctions. 

Let us now consider the Green's functions in the regulated field theory at finite N, 
Gn{xu ...,xr;yi,...,yr) = C {^\x^) $(yi) ■ ■ • , (7.2) 

where the parameter = (^tv > is a multiplicative wavefunction renormalization. Then 
the regularized partition function ()6.15|) and all Green's functions ()7.2|) at finite are 
invariant under the scale transformations 

X I — > A X , 
B ^ X-^B, 

e I — > x^e, 

2 \ -2 2 

/i I > X jjL , 

C ^ A"-iC, (7.3) 

where we have redefined coordinates using (j7.ip so that the scaling parameter s = 1 in 
fl6.7|) - ()6.1U|) . This property is elementary to check by using (17.11) and compensating the 
transformations ()7.3|) by the change of matrix integration variables M — >■ A M, M"^ — > 
AM^ The Jacobian A^^ of this latter transformation is not important here, as it cancels 
out in the Green's functions and leads to a finite, irrelevant shift of the free energy. 

This scaling property is important for the physical interpretation of the large limit 
of the matrix model. It implies that we are free to make the transformation ()7.3|) with 
A = A^'' for arbitrary real u, before taking the limit N oo. The result is independent 
of the exponent u. However, even though the mathematical expressions for the Green's 
functions are unchanged, the physical interpretation of the large limit does depend on 
u. Choosing the scaling x — ^ N'^ x and sending N oo for u > means that we blow 
up small distances and "zoom" into space. On the other hand, for u < the opposite 
phenomenon occurs, we "zoom" out so that short distance structures become invisible. 
The first scaling limit is an infrared limit while the second one is an ultraviolet limit. To 
be more precise, let us put the system on a hypercubic lattice of side R and spacing a, 
so that there are A^ = (-R/a)^" lattice points. In conventional quantum field theory, one 
usually performs the ultraviolet limit, which means sending a ^ keeping R fixed (and 
adjusting all parameters so that the limit makes mathematical sense). In the infrared (or 
thermodynamic) limit, one sends i? — >■ oo keeping a fixed (so that the finiteness of the 
spatial volume no longer matters). 
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This suggests that the parameter A introduced in ()3.24j) has different interpretations 
dependent on the choice of exponent v. In the two cases described above we may identify 
A oc with a— >^OorAcxl/a with R —>■ oo, respectively. There can also be many 
other interpolating identifications as well. The definition fl3.24|) is simply the statement 
that there is a length scale in the field theory, the renormalized magnetic length. Any field 
theoretic interpretation of this length is possible, and this point of view is consistent with 
the duality described in Section 12.31 This is one of the various unconventional properties 
possessed by the present class of quantum field theories. In the ensuing subsections we 
shall make these statements precise and explicit. 

For this, we will introduce the two basic cutoffs above more precisely. The infrared 
cutoff R is like the diameter of space and is defined by 

i?2 = 47riV' e , (7.4) 

with A^' 7^ A^. The significance of this length scale is well-known from the physics of the 
quantum Hall effect. To understand its interpretation better, consider the i = 1 Landau 
eigenf unctions which from ()3.8|) and ()3.9p are given explicitly by 

^m— 1 

01 Ux) = , e -l^l'/^e . (7.5) 

^ v/tt^ (m- 1)! ^ ^ 



The one-particle "location" |0i,m(a^)P thereby has a sharp peak at \z\ oc ^J{rn — V) 9, and 
so restricting the quantum numbers m < A^' is equivalent to the bound \z\ < R above. 
This argument also works for i > 1. In other words, the infrared cutoff ()7.4|) amounts to 
constraining the system quantum numbers to lie below a finite number A^' within each 
Landau level, as in the quantum Hall effect. This cutoff was ignored previously and 
implicitly sent to infinity. 

The other cutoff is our standard ultraviolet cutoff which is the lattice spacing a defined 

by 

" = A^ = AT • 

For clarity, in the remainder of this section we shall work mostly in two dimensions n = 1 
and only with the finite A^ two-point function which is given by 

TV oo 

Gn{x, 1/) = - X] ^^(^^) 'PrnA^) <Pi,m{y) ■ (7.7) 

e=i m=l 

Here ei = AB{t — |) -|- //^ and W]y{^) is a function which converges to the rescaled loop 
function W{C,/A'^) of the matrix model in the limit A^ — ^ oo. According to ()7.3p we 
may set C = 1 in two dimensions. Note that ee ^ + iJ? at large A^. 



7.2 Infrared Limit 

The particular scaling we studied earlier, whereby a finite ultraviolet cutoff is kept, will 
be refered to as the "infrared limit" of the noncommutative field theory. Explicitly, it is 
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defined by first taking R ^ oo, then N ^ oo such that a is finite. Looking at the non- 
interacting Green's function ()5.13|) . the asymptotic behaviour Kq{z) ~ — In^r for z —>■ 
imphes that this limit corresponds to "zooming" into space and the breaking of translation 
invariance by the magnetic field becomes invisible. In terms of the scale transformations 
of the previous subsection, it corresponds to setting u = —1 and rescaling the couplings 
according to (j7.3|) with A = I/a/^tt^. The quantity (j3.24j) is thereby interpreted as an 
ultraviolet cutoff and the Green's function is computed in the limit 

Gir{x,y)= lim GN{x,y) (7.8) 



N- 



=Af/47rA2 



keeping A finite. The result is given by ()5.21|) . with — > ^ in the limit. Everything 
is known about this limit. The quantum field theory is quasi-free, i.e. everything can be 
calculated solely from the explicit knowledge of the two-point Green's function. There is 
no new length scale introduced between a — > and oo. 



7.3 Ultraviolet Limit 

Combining the result above with the duality described in Section IT^ immediatelv implies 
the existence of another scaling in the noncommutative field theory in which the physics is 
completely different. We will refer to this dual scaling limit as the "ultraviolet limit". In 
fact, it is straightforward to establish the duality property of the finite N Green's function 
()7.7|) for fixed parameters g and 6. This follows directly from the generating functional 
f|3.8|) . which can be easily checked to possess the symmetries Vs,t{—x) = V-s-t{x) and 
'Ps,t{p) = 2'n'9Vis,it{0p)- It follows that the Landau eigenfunctions obey 

4>,,Up) = -27ri^+"^^</)„/ep) , 

-x) = (-l)^+'"0,,^(x) , (7.9) 



and hence that the regulated propagator has the duahty 

G^{p,q) = i27r9y GM{9p,eq) (7.10) 
under Fourier transformation. 

This property implies that the limit 

1 ~ f p q 



G^Mq)= \un^—G^ (^'^ 



(7.11) 



=Ar/47rA2 



keeping A finite is also well-defined. By using ()5.21|) and substituting x = — p/47rA^ it 
can be written as 

1 / p q 



Guv{p,q) = 



47rA2 ' 47rA2 



^ J d\ W^b((47rAa;)2) e-^"-(P-'') , (7.12) 



A2 

\x\<l/y^A 
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and it thereby gives an ultra-local Green's function in the ultraviolet scaling limit/ 

GUx,y) = -^W^b((47rAx)2) ^^'^(x-y) (7.13) 
with |x| < I/a/tFA. In particular, in the free case it reduces to 

It is interesting to note that this scaling limit could also have been defined directly as 

GUx,y)= lim N^GN{Nx,Ny) . (7.15) 

The limit in ()7.15|) is technically rather difficult to carry out explicitly. But the duality 
symmetry of the quantum field theory enables a quick and easy derivation from the 
corresponding results of the infrared limit above. 

This scaling corresponds to first taking the limit a — 0, then N' ^ oo such that R 
is finite. In terms of the scale transformations of Section 17.11 it corresponds to setting 
= 1 and rescaling the couplings according to ()7.3|1 with A = AttO. Note that this 
means now ^ — >■ and A is interpreted as a finite infrared cutoff. From (|5.13|) and the 
asymptotic behaviour Kq{z) ~ for — > oo it follows that this limit corresponds 



to "zooming" out of space, i.e. probing large distances of order N/ K = y/AnN 9, and all 
that is observable is the spatial dependence due to the magnetic field whereas the center 
of mass dependence becomes ultra-local. Because of the UV/IR duality, this ultraviolet 
limit is given by the same 't Hooft limit of the matrix model, and all that happens is 
that the long and short distance cutoffs of the theory are interchanged. Everything about 
this limit is also understood, because its properties are just dual to those of the previous 
subsection. In particular, the quantum field theory is again quasi-free, but with two-point 
Green's function given by the ultra- local form ()7.13p . and there is no new length scale 
introduced between and R oo. 

While this argument for going from the infrared to the ultraviolet limit is mathemat- 
ically very simple, it is somewhat paradoxical. In the infrared limit we take 6 oo and 
then, by a simple rescaling of the spatial arguments leaving the matrix integral unchanged, 
we arrive at a limit whereby 6' — > 0. The resolution of this paradox comes from using 
the invariance of the matrix integral under M — >■ N'^ M, M'^ — > N'^ M'^ before the limit 

— > cxD is taken. Thus if we consider just the matrix model, then from (j7.3|) it follows 
that the only relevant parameters are 9 and 9g. This suggests a different physical in- 
terpretation of the infrared limit above which makes the existence of the ultraviolet limit 
somewhat less paradoxical. The crucial computation for the infrared limit was the sum 



^Note that neither the ultraviolet nor the infrared limits respeet magnetic translations. The reason was 
explained at the beginning of Section [2. 21 Magnetic translation invariance implies translation invariance 
in both position and momentum space (up to phase factors). However, the limits 6* — > cx) and 6* ^ are 
singular. In the latter case the translational symmetry is broken in position space (and the propagator is 
ultra-local there), while in the former limit it is broken in momentum space. The breaking of translational 
invariance also occurs to all orders in perturbation theory of standard noncommutative <i>'*-theory with 
a finite cutoff 0. 
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over Landau levels ()5.19p for large ^ (of order A^). However, following the derivation of 
Appendix G through, it is clear from ()7.7p that the same result can be obtained also in 
some limit keeping the noncommutativity parameter 9 finite and 



lim 

Af— >oo 



G 



N 



X 



(7.16) 



6»=l/47rA2 



The results obtained are exactly as before, only now A has a different interpretation. It 
is simply the inverse of the length scale set by the noncommutativity parameter, and it 
need not be taken to infinity any longer. In a similar vein, we may write 



lim 



NG 



N 



Nx, VNy 



= l/47rA2 



(7.17) 



This makes clear the interpretations above as zooming in versus out of space for the in- 
frared versus ultraviolet limits. It is natural from the point of view of the noncommutative 
duality. 

This discussion motivates the search for alternative scaling limits in between the other 
two that we have analysed, in which the Green's functions have a non-trivial structure 
on both short and large distance scales, and whereby the magnetic translation symmetry 
is preserved. The simplest one uses both infrared and ultraviolet cutoffs ()7.4|1 and ()7.fij) . 
and corresponds to a scaling exponent 1^ = with A = 1. It comes from sending a 
and — *• cx) such that 6 is finite, i.e. N' = N ^ oo with 6 finite. This can be regarded as 
a "true" quantum field theory limit, in which both ultraviolet and infrared cutoffs have 
been removed. It is defined by 



Go{x,y) 



lim 

N^cx 



GN{x,y) 



9=l/47rA2 



(7.18) 



and it thereby corresponds to keeping all length scales in between the infrared and ul- 
traviolet regimes (no "zooming" in or out). We believe that such intermediate scaling 
limits will only lead to an interacting quantum field theory if some non-planar limit of 
the matrix model is used. Such an analysis is beyond the scope of the present work. 



7.4 Lifting the Degeneracy of Landau Levels 

The alternative scaling limits of the noncommutative quantum field theory that we have 
proposed in this section, and in particular the introduction of the intrinsic infrared cutoff 
()7.4|) . have the effect of essentially removing the degeneracies from each Landau level. As 
is known from the theory of the quantum Hall effect, an alternative way to accomplish this 
is to add a confining electric potential to the Landau Hamiltonian, as we did in (j2.1|) . Thus 
this model has the potential of also providing an exactly solvable interacting quantum field 
theory. There are several ways that one may proceed in analysing this perturbation of the 
original field theory. For instance, one may regard a in ()2.H) as a small parameter and 
simply examine the perturbative corrections in a to the results above. This is tractable, 
as we have illustrated how to completely solve the a = model. Alternatively, one could 
pursue a saddle-point analysis of the eigenvalue representation ()4.6|) in the large limit. 
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The effective eigenvalue action describes interacting particles in a common self-interaction 
potential V but with an electric field dependent on each individual particle. Particle i 
feels an electric field i/6. Thus the equilibrium configuration of the particles occurs when 
they are well-ordered. Permutation symmetry is then broken, and rather than A^! identical 
saddle points, there is just a single one. This is true for very large a, i.e. for equilibrium 
positions hi such that a\hi — h^l ^ 0- For small a the eigenvalues hi with i ^ 1 
accumulate close to the origin and the logarithmic Vandermonde type repulsions at short 
distances become important. Finally, one can generalize the auxilliary Hermitian matrix 
integral representation of Appendix C to this model, and analyse it as described there. 
The solution of this generalized noncommutative field theory represents an interesting 
challenge in the search for non-trivial, exactly solvable interacting models. 

Another, more difficult way to lift the degeneracy of the Landau levels is to perturb 
the field theory away from the special point B = 9~^. This of course destroys its GL(oo) 
symmetry and hence its exact solvability. From ()3.18p we see that perturbing the action 
to B ^ 9^^ is like adding hopping terms, giving a model defined on a two-dimensional 
lattice with somewhat unusual shifting actions. One way to proceed in solving this 
lattice model would be to try to map the terms involving the shift matrix ()3.20j) onto 
the generators of an auxilliary Heisenberg algebra represented by infinite-dimensional 
matrices. The total kinetic energy term would then assume a more natural form. It 
may be that these additional hopping terms produce relevant deformations of the trivial 
Wilsonian renormalization group fixed point which we have found. A similar sort of 
perturbation expansion has been studied in P OHj. 

Coupling scalar field theories to fermionic field theories with enough supersymmetry is 
known to be a way to rid the noncommutative quantum field theory of UV/IR mixing [221 
1^ . It may be that the fermionic analog of the quantum field theory with action (j2.im) 
contributes power-like divergences which cancel those of the bosonic theory, and therefore 
that appropriate supersymmetric extensions yield a renormalizable interacting model. 
Such supersymmetric field theories are the topic of the next section. 



8 Supersymmetric Extension 

In this final section we will formulate a supersymmetric extension of the scalar quantum 
field theory defined by the action ()2.10j) . We will only present its explicit construction 
in detail, and leave aside its exact solution for future work. The particularly interesting 
aspect of this formulation is that non-trivial supersymmetric interactions can only be 
formulated using the noncommutative star-product. This fact is even apparent, as we 
will see, at the level of the free supersymmetric action, whereby the supersymmetry 
transformations are intrinsically operator- valued and hence are parametrized by elements 
of a noncommutative algebra. We will begin by writing down the analog of (j2.H) with 
CT = 1, a = for relativistic fermions, and then proceed to describe the supersymmetric 
combination of the two field theories. Throughout this section we work in two spacetime 
dimensions and at the special point B = 6^^. 
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8.1 Fermionic Models 



Consider the noncommutative field theory of a massive two-component Dirac fermion field 
"^{x) in two dimensions under the influence of a constant background magnetic fleld. The 
action is 

Si = J d^x ^^(x) {a' Di + fi^) ^(x) + K (^^ * ^) (x) , (8.1) 

where a\ i = 1,2 are the usual Pauli spin matrices. The Dirac operator can be diago- 
nalized in the complex coordinates ()3.1|) by expressing it in terms of the ladder operators 

as 

.'D. = 4r. 'A . (8.2) 



In terms of the Landau eigenfunctions, the orthonormal solutions of the Dirac equation 
are then given by 

= ( tr'\l.) (8.3) 

with 



^/2 \±i 0f,m(x) 



^^D,^,=^^ = ±2iy^^J^, (8.4) 

where we have used (j3.17|) . Here m = 1, 2, . . ., while i = 2,3, . . .. 

The case of £ = 1 should be treated separately. The ± eigenfunctions ()8.3|) are 
interchanged under the action of the two-dimensional chirality operator, which is the 
diagonal Pauli matrix in the basis of Dirac matrices that we have chosen. In contrast, 
zero modes of the Dirac equation are all chiral, as only one independent solution for each 
m exists. They are given by 

^l,m(x) =( .] (8.5) 



with 

a^D,^i,^ = 0. (8.6) 

The existence of chiral zero modes of the Dirac operator follows from the index theorem, 
according to which the number of zero modes of negative chirality minus the number of 
zero modes of positive chirality is equal to the total magnetic flux divided by 2tt. In the 
present case the total flux is inflnite, as is the number of chiral zero modes. 

We can now expand the fermion fleld of (|8.ip in the complete basis of Landau wave- 
functions to get 

oo / oo \ 

^{x) = {Aney/' J2[Y.J2 ^^^^^ ^^'-(^) + /™ ^o,n.(^) , (8.7) 

m=l \ e=2 s=± / 
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where F± = {F±^rne)m>i ,e>2 are complex Grassmann-valued matrices, and / = (/m)m>i is 
a complex Grassmann-valued vector. Note that F± are rectangular matrices, in contrast 
to the bosonic case. In terms of them, the action ()8.1|) can be expressed as an infinite, 
complex fermionic matrix model coupled to an infinite, complex fermionic vector model 
by using the star-product identities 

(^l,m)^*V'£'',m' = S"'' 6u> (f)m,m' , 

1 



s 

m 



which follow easily from ()3.12|) . The action thereby becomes 

= 5^ Tr^(F,/C^Ft)+5^v^;Xf Tr^(FtF,)+v^^f/t/ 




.9) 



where 

]C% = ± i Vl67r(£-1) 6u' . (8.10) 



8.2 Supersymmetry Transformations 

We will now combine the fermionic model of the previous subsection with the bosonic 
model of Section 12.11 and define a supersymmetric field theory on a noncommutative 
phase space. At a first glance, it appears impossible to have supersymmetry in such a 
field theory, because bosons and fermions interact with an external magnetic field in a 
different way. Indeed, the spectrum of the free bosonic Hamiltonian, 2B{2i — 1), is offset 
by the quantity 2B with respect to the spectrum of the square of the free Dirac operator, 
4B{i — 1). However, this offset can be compensated by a shift in the boson mass. To see 
how this works, let us compute the partition functions for free scalar fields of mass /x and 
free fermion fields of mass /if. Up to irrelevant constants, they are given by 



'(0) 



£,m=l 
oo 

n 

m=l 



/if n + - 



(8.11) 



If we now take fi^ = fi^ — 2B, then the functions in (|8.1ip are almost the same, up to a 
minor mismatch in their zero mode contributions. 



It is straightforward to see that the free action for the complex scalar field and the 
Dirac fermion field, with masses related as above, is indeed supersymmetric. Using ()3.16j] 
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at B9 = 1 and ()8.2|) . the action can be rewritten in terms of ladder operators as 



'susy 



j d'x ^\x) (A+A + /i^) $(x) + ■^\x) ^ ^ ^(x) (8.12) 

where A = 2 ^/B a. The action (|8.12p is invariant under the infinitesimal supersymmetry 
transformations 



A<|.t = 



fi^ 

A 

/if 



where the parameters e and of the supersymmetry transformations are arbitrary Grass- 
mann odd functions of the b and operators in ()H.2j) . Here we use the convention that the 
action of a differential operator from the right is consistent with integration by parts, i.e. 

/a^bt/, /at^b/. (8.14) 

It follows that the supersymmetry transformations in this case contain infinitely many 
parameters. This resembles local supersymmetry somewhat, except that here the param- 
eters of the transformation are arbitrary functions of differential operators rather than 
arbitrary functions of the spacetime coordinates. A closer analogy is the type of supersym- 
metry that arises in zero-dimensional supersymmetric matrix models [Ml^ISni; in which 
the parameter of the transformation is an arbitrary matrix. In fact, the parameters of the 
supersymmetry transformation in the present case become matrices after expanding in 
the basis of Landau eigenfunctions. For this, we note that the kinetic energy operator for 
the boson field in (|8.12j) can be written in matrix form as )C~^)C~ , where /C^ are the ma- 
trices defined in fl8.1()|l . The action ()8.12|1 can thereby be written as the supersymmetric 
matrix-vector model 



C(0) _ ry 



J2Fs (jC' + Vi^/if) F| + (/C+/C- + 4n9n^f) 



+ fi, f f + AttO [fi^, - 20-^) (3^ (3 , (8.15) 

where for convenience we have separated out the zero-mode part = ^im of the scalar 
field, and in ()8.15|) it is understood that M = {Mim)e>2,m>i is a rectangular matrix. The 
action ()8.15p is invariant under the infinitesimal supersymmetry transformations 



6,F} = - (^IC-' + Vi^/if) el , 

s=± 

5,Mt = 5^F,et, (8.16) 



s=± 
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where e± = {e±/ei)e^i'>2 is a pair of infinite dimensional matrices. As expected, the 
supersymmetry transformations do not involve zero modes. 

Since /C"*" = {JC~y, we have {S^FsY ^ ^eF}, and hence and should be regarded 
as independent variables in the functional integral defining the corresponding quantum 
field theory. The supersymmetry transformation (j8.16|) is then a legitimate change of 
variables, and it leads to the supersymmetric Ward identities in the usual way, as long as 
the integration measure is superinvariant. The condition for this is 



The absence of zero modes in ()8.16|) and the constraint ()8.17|) eliminate only a finite 
number of degrees of freedom. 

8.3 Supersymmetric Interactions 

Thus far we have considered a non-interacting supersymmetric quantum field theory for 
which noncommutativity played no role. Noncommutativity becomes important if we 
try to answer the question of whether or not there exist interactions that preserve the 
supersymmetry. We believe that the affirmative answer to this question can be given only 
within the framework of noncommutative field theory. It seems impossible to construct 
local interactions that are invariant under the supersymmetry transformations (j8.13|) . 
though we cannot prove this fact rigorously. However, "star-local" interacting quantum 
field theories which possess the infinite-parameter supersymmetry (j8.13|) do exist, as we 
now proceed to demonstrate. 

The reason why noncommutativity aids the construction of Lagrangians which are 
invariant under the supersymmetry transformations of the previous subsection is that the 
ladder operators (|3.2|) obey remarkable Leibniz- type rules with respect to the star-product. 



These equalities can be proven by expanding the fields / and /' in the Landau basis 
and using ()3.12p . Alternatively, the star-product projector relation ()3.12|) for the Landau 
wavefunctions can be derived from the Leibniz rules (|8.18|) . which can be independently 
checked by a direct computation in the Fourier basis by using (j2.5|) . Formally, the relations 
(|8.18|) simply refiect the fact the algebra of functions on M^, equipped with the star- 
product, generates a bimodule for the harmonic oscillator algebras (j3.3p realized by the 
commuting a and b operators. This was already implicit in the definitions ()8.14|) . 

As usual, the construction of supersymmetric Lagrangians is facilitated by the intro- 
duction of auxilliary fields JF. With them, we postulate the supersymmetry transforma- 
tions 




(8.17) 



a(/^/') 

b(/^/') 
(/ ^ /') 
(b/)^/' 
(bV)^/' 



(a/)^/', 
(aV)^/', 
/^(bf) , 
/^(bV) , 

/ ^ (a^ n , 

/^(a/) . 



(8.18) 
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(8.19) 



where C is an arbitrary parameter of mass dimension 1. It can be set to unity, but then 
the left and right components e and et will have different scaling dimensions. 

There are two star-quadratic invariants of the supersymmetry transformations ()8.19j) 
which contain at most two derivatives of the fields. They are given by 



vl/t 



A 

-At 

^t ^ ^ + i jrt ^ $ + i $t ^ jr . 



^ + $t*AtA$ + J^t^jr 



(8.20) 



The superinvariance of the expressions ()8.20j) depends crucially on the bimodule prop- 
erties ()8.18p of the ladder operators. We can now proceed to construct supersymmetric 
Lagrangians by taking star-products of L\ and £o- The most general renormalizable 
Lagrangian, i.e. the one with no couplings of negative dimension, is then given by 



(8.211 



The auxilliary fields can be eliminated from (j8.2H) by using the equations of motion. Then 
the second term gives masses to the boson and fermion fields, while the third term induces 
various interactions. These interactions are star- local, but non-polynomial in $ and 
Explicitly, one finds 



5-. 



susy 



d^X 



-At 



A 



$t(x) (AtA + /i^) $(x) + ^t( 
-f 2gii\ $t ^ $ ^ $t ^ $(x) 2giii $t ^ $ ^ vpt ^ ^{x 



^(x) 

^v[/t^\I>^vl>t^vl>(a;) 

2 ^ ^ 



(8.22) 



The action (j8.22|) possesses the usual star-local GL{oo) symmetry as before and hence 
yields a potentially solvable supersymmetric quantum field theory. It would be interesting 
to seek a superspace formulation of this model, which may aid in finding its exact solution. 
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However, it is not clear what the concept of a superspace could mean in the present con- 
text. The supersymmetry here acts as a fermionic rotation on the Landau levels, which 
has nothing in common with ordinary supersymmetry. In particular, the commutator of 
supercharges is not a spacetime translation . The definition of noncommutative 

superspaces has been addressed recently within various different contexts in |SZj-[nn]- In 
the present case, supersymmetric rotations are the super-analogs of the bosonic GL{oo) 
symmetry, so it is tempting to speculate that together they generate an infinite dimen- 
sional GL[oo\oo) supergroup. Geometrically, this noncommutative supersymmetry would 
then correspond to the superspace generalization of area-preserving diffeomorphisms. 
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Appendix A Mean Field Analysis of Commutative 

Complex <l>^-Theory 

In this appendix we point out that the models studied in this paper also have a concrete 
physical motivation as a novel kind of mean field theory for conventional (commutative) 
complex $^-theory. For this, we consider the action 5*0 + ^int for ordinary two-dimensional 
complex $^-theory written in momentum space, with free part 

5-0 = ^ d^k {e + /i') ^\k) l>(fc) , (A.l) 

and the interaction 

Sint = y J d^k d^k' d'g d\' S'-^\k' - k + q' - q) ^\k') ^k) ¥{q') $(g) (A.2) 

describing all possible momentum conserving processes whereby particles with incoming 
momenta k and q are scattered to new momenta k' and q'. Standard mean field theory 
can be applied to this theory by truncating the interaction to retain only the Hartree and 
Fock terms. In the former interactions the individual particle momenta are conserved. 
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k = k' and q = q', while in the latter terms they are exchanged, k = q' and q = k'. In 
the present case these two types of terms are actually the same, and their sum yields the 
Hartree-Fock interaction 

Shf = 90 j d^k d'q ^^(k) ^k) <|t(g) <|(g) . (A.3) 

The quantum field theory with action 5*0 + ^hf can be solved exactly, and its exact 
solution is identical to mean field theory for commutative complex $'^-theory. Up to now 
our discussion could have been easily extended to arbitrary spacetime dimensionality, but 
in two dimensions there are special mixed interactions which are of Hartree type in one 
component of the momenta and of Fock type in the other component, k' = {k'i^k'2) = 
{ki, ^2) and q' = {q[, gg) = ^2), or vice versa j^Zl. The sum of all of these mixed terms 
yields the interaction 

^mixcd = 90 j d^k d\ ^\ki, q2) $(A;i, k2) ^\qi, ^2) ^2) • (A.4) 

We now claim that the model with action 5*0 + ^mixed is equivalent to the one studied 
in this paper. For this, we introduce a regularization by restricting the momenta to 
k = {ki,k2) = ^{i,m) with i,m = 0, ± 1, ± 2, . . . , ± i?/a = (iV- l)/2, where i? < 00 is 
a large but finite radius for the size of spacetime serving as an infrared cutoff, and a > 
is a small but finite lattice spacing serving as an ultraviolet cutoff. We can then identify 
the Fourier modes of the fields with the elements of an x iV matrix M = (M^m) as 

/^27rV /27r£ 27rm\ 
M,^ = — $ — , -— . (A.5) 



J ^ \ R ' RJ 
By introducing a diagonal matrix E with elements 

Etn. = e hm , (A.6) 

we may thereby write the regulated action as a complex matrix model 

^"0 + ^mixcd = Tt (^(jM^ EM + dMEM^ M + | M^MM^M^ (A.7) 

with a = a = 1, g = 2go/R'^, and Tr the usual N x N matrix trace as before. 

The matrix integral that we evaluate in Section 14.11 provides an exact and explicit 
formula for the free energy of this model. In Sections ID and we find the complete 
solution of the model in the limit R ^ 00 for the special maximally anisotropic case 
a = 0. In section [7| we examine other possible ways of removing the cutoffs in the 
anisotropic model. The physical interpretation of the model presented here gives a strong 
motivation for extending the results of this paper to the more complicated isotropic case 
where a = a. 



Appendix B UV/IR Duality versus T-Duality 

Combining the translation and symplectomorphism invariances of the action (j2.1Up with 
its Morita-type duality symmetry found in [T2] lends further support to the observation 
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of ^21 that the quantum field theory defined by ()2.10|) may be regarded as a discrete 
noncommutative Z12 gauge theory when B = 6~^. For this, we introduce a constant, 
background metric H = {H^y) on R^" and use it to rewrite the action ()2.1U|) in a covariant 
form. We regard the action as a functional of the dynamical field $, the background 
fields H and -B, and the coupling parameters g and 9. The action then has a duality 
under Fourier transformation [T^ . 



S^[^]H,B,g,d] = det{B/2n 



^;H,B,~9,e 



(B.l) 



where the dual parameters in momentum space are given by 



H 
B 

9 
6 



B-^HB-^ 
9 



det(27r^ 



(B.2) 



Note that, in contrast to the transformations of Section 12.31 here we do not rescale the 
arguments of the fields. 

There is a novel gauge Morita equivalence interpretation of this noncommutative du- 
ality, wherein we may heuristically regard the model ()2.10|) as a noncommutative gauge 
theory defined along 2n finite discrete directions, i.e. on a two-sheeted manifold M^" x Z2. 
Noncommutative Yang-Mills theory defined on a 2n-dimensional torus is manifestly in- 
variant under the standard SO{2n,2n,Z) open string T-duality transformations ^ 

e = {Ae + By {OQ-N)-^ , 

~9ym = 9ym det(^Q-N) , (B.3) 

where the superscript T denotes transposition, gyu is the Yang-Mills coupling constant, 
and in this equation 6 is the dimensionless noncommutativity parameter. The 2n x 2n 
symmetric integral matrix N is proportional to the rank of the gauge theory, while the 
antisymmetric integral matrix Q is determined by the magnetic fluxes of the gauge bundle 
around the various two-cycles of the torus. The symmetric and antisymmetric integral 
matrices A and B, respectively, are chosen to solve the generalized Diophantine equation 

AN + BQ = l2„. (B.4) 



Going back to the scalar field theory we may, in the usual Connes-Lott type inter- 
pretation of quartic scalar field theory [SH' , regard $ as the off-diagonal components of a 
super connect ion, and the field theory ()2.10|) as induced by the limit in which the diagonal 
gauge degrees of freedom are frozen out in the usual noncommutative Yang- Mills action 
for the superconnection. The noncommutative scalar field theory is then formally the zero 
rank limit of a noncommutative gauge theory. All formulas above still make mathematical 
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sense in this limit and it is heuristically the choice that should be made for a discrete 
gauge theory. It means that one should set 

N = . (B.5) 

The equation ()B.4|) is then solved by 

A = 0, Q = -B=(^° "p^^^ln. (B.6) 

Then the transformation rules ()B.3|) coincide exactly with ()B.2|) at the special point 
B = 9~^. This holds with the identification (yfyM ~ ((^tt)" f?)"^. For this "discrete" 
truncation of the toroidal compactification, the choice of flux matrix Q in flB.6|) is very 
natural, given the charges carried by the scalar fields $ (see (I2.33p ). This gauge theoretic 
interpretation also indicates why the duality property is uniquely possessed by a quartic 
interaction of the scalar field, which has other special features as well which are displayed 
in the main text. It would be interesting to pursue this gauge theoretic interpretation 
further using the usual techniques of noncommutative geometry, as it may help elucidate 
further features of the model and its generalizations. 



Appendix C Alternative Derivation of the Loop 

Equations 

In this appendix we will present another derivation of the Schwinger-Dyson equations 
fl4.25|) for the matrix model ()3.26|) with potential ()2.9|) . This derivation has the advantage 
of being somewhat more amenable to generalization to other noncommutative complex <l>^ 
field theories, such as that associated with the matrix model ()4.1|) . It also gives another 
physical interpretation to the loop function (|4.23|) in terms of the spectral distribution of 
an auxilliary Hermitian one-matrix model. 

We introduce an auxilliary N xN Hermitian matrix X = (Xim) through the Hubbard- 
Stratonovich transformation 

where DX = Yli dX^^ nm<£ ^ ReX^m d ImX£m/27r. After a rescaling of X and M by the 
matrix dimension N, the partition function ()3.26|) may then be written in the form 

Z^{E) =JdM DMt y DX e -i ^'-MHmE- ^ixm)M ^ ^^ 2) 

where for convenience we have changed variables from N x N complex matrices M = 
(Mim) to vectors M G through the rule 

M = (Mn, M12, . . . , Mi^, M21, . . . , M^^y . (C.3) 



57 



The Gaussian integration over M and M^' can now be carried out in this manner, giving 
the effective Hermitian one-matrix integral with a Penner type interaction potential, 



Zm{E)= j DX e-iTVX2-Ti-^,ln(10i?T+iX®l) 



where Tr is the matrix trace in the adjoint representation. 



(C.4) 



The matrix integral ()C.4|) is invariant under unitary transformations of both X and 
and can thus be written as an eigenvalue model. For this, we diagonalize X = 
U dia.g{fii, . . . , ^in) U\ with /i£ G M and U G U{N), compute the Jacobian of the in- 
tegration measure to get DX = [dU] Yle^f^e ^-^Nlf]'^, and thereby write ()C.4|) as 

Zn{E) = n / d/i, e-i'^^' n ] . ■ (C.5) 



1=1 



The saddle point equation for the large N limit of the eigenvalue integral ()C5|) is 

We now multiply the expression ^C^\ through by {fig — z)^^ for z G C, sum over all 
i = 1, . . . ,N, and expand the terms to rewrite them in terms of the resolvent function 

1 ^ 1 

By rewriting the sums over the eigenvalues of the external field E as integrals over 
the spectral density ()4.24|) . and dropping the term T,'{z)/N a,t N = oo, we may in this 
way write ()C.6|) as the large N loop equation 



a2 

\{l + z S(^)) + T?{z) + jdX (S(^) - S(- i A)) = . (C.8) 

The loop equation ()C.8|) is essentially the same as ()4.25|) . Since the resolvent function 
flC.7|) has the asymptotic behaviour T,{z) ^ — ^ + O(^) for z oo, it then follows that 
it may be identified with the function ()4.23|) through 

J:{z) = -iW{-iz) . (C.9) 



This coincidence is not surprising, as both functions W and S are the generating 
functions of correlators of derivatives of M'^M acting in the complex matrix model. Indeed, 
a straightforward calculation using the Hubbard-Stratonovich transformation above shows 

fe, wiX,) . 1 ((MtM)„,)^ = <J 1 ± ^-1—^ , (C.IO) 
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where the second expectation value is taken with respect to the Hubbard- Stratonovich 
representation ()C.4|) . The usefulness of this approach, however, is that the eigenvalue 
model ()C.5|1 may be straightforwardly generalized to include also the external field E 
that appears in the matrix model ()4.H) . and its loop equations may be analysed by the 
technique presented here. The exact solvability of the matrix model ()4.1|) may thereby be 
attributed to the U{N) symmetry of the NxN matrix model of the Hubbard- Stratonovich 
field X. This then has the potential of presenting an exact, non-perturbative solution of 
the generalized noncommutative quantum field theory with action and interaction 
potential given by ()2.9p . Furthermore, this technique has the advantage of straightfor- 
wardly supplying an explicit formula for the generating functional ()2.13|) of the connected 
Green's functions of the model. 



Appendix D Explicit Solution of the Master Equa- 
tion at Large N 

To solve the non- linear integral equation (j4.25|) . we will reduce it to a standard Riemann- 
Hilbert problem of the type that appears for ordinary Hermitian one-matrix models. For 
this, we introduce the resolvent function of the external field E which is defined as the 
Hilbert transform 

= /dA (D.l) 

ai 

of the spectral density ()4.24|) . and the function 

= ?dA . (D.2) 

J z — X 

Both of these functions are analytic everywhere in the complex z-plane, except on the 
support [01,02] of the spectral distribution where they each have a branch cut. Let 

cu±(A) = i (^cj(A+ iO) ±cj(A- iO)) , AG [01,02] (D.3) 

be the continuous and singular parts of the function uj{z) across its branch cut. We simi- 
larly define fi± (A) for A G [01,02]. Note that co'_(A) = i7rp(A) andfi_(A)= inW^X) p{\). 



From the definitions ()D.1|) and ()D.2|) it follows that the singular parts of these functions 
are related through 

n^X) = W{X) uj4X) , VA G M . (D.4) 

Furthermore, from the large Schwinger-Dyson equations ()4.25p it follows that the con- 
tinuous part of ()D.2|) obeys 



~gn+{X) = ~gW\X)-(^X-~gu;+{X)')w{X)~l , Ag[oi,02]. (D.5) 
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Assuming g ^ 0, this suggests the analytic ansatz 

n(z) =W\z) - ( - -uj(zy\ W(z) , yzeC. (D.6) 

/ 9 

Substituting ()D.6|) into ()D.4|) and ()D.5|) then imphes the respective restrictions 

W^{\) (^^-uj+{\)-2W+{X)^ = 0, VAgM, (D.7) 

W.{X)(W-i\)+U-{\)^ = 0, \e[ai,a2] (D.8) 

on the real-valued analytic function W{z) in the complex z-plane. In deriving ()D.7|) from 
(ID.6|) we have assumed that W{z) has no poles. 



In addition to ()D.7|) and ()D.8|) . a third restriction on the function W{z) is imposed by 
the boundary conditions required to solve the Schwinger- Dyson equations. From ()4.25j) 
it follows that there are two branches of solution, one with the asymptotic behaviour 
W{z) ~ I for z oo, and the other with W{z) ^ — ^ for z — oo. We will take the 
latter boundary condition, as it is the one which matches that of the perturbative solution 
(I4.27P of the complex matrix model. This branch reduces continuously to the Gaussian 
solution W^'^\X) = — J at g = 0.^ Thus we will in addition require 

W{z)c^-- + o(\] asz^oo. (D.9) 



z \ z^ 



We now solve ()D.7|) - ()D.9|) by making a one-cut ansatz, which will lead to a rational 
parametrization of the solution. For this, we assume that W-{X) is non- vanishing on a 
single connected interval [61,62] in the complex A-plane. To satisfy ()D.8|) we must then 
have 

[ai, 02] n [61,62] =0 . (D.IO) 

Since the resolvent uj{z) is analytic everywhere away from its branch cut on [ai, 02], from 
(lmn|l it follows that u;_(A) = for A G [61, 62]. From (lD?7j) we then have 

W4X) = ^-uj{X)-^ , A €[61,62]. (D.ll) 



We can now turn ()D.11|) into a standard Riemann-Hilbert equation 

A €[61,62]. (D.12) 



WiX) \ _ t-^(A) 



v/(A-6i)(A-62) / 2i ^(62-A)(A-6i) 



^The former asymptotic behaviour is the pertinent boundary condition to use in solving the corre- 
sponding Hermitian matrix model in an external field, obtained by the formal substitution X = Af^M 
in the matrix integral (|3.26|) . along with the appropriate change of integration domain. This Gaussian 
Hermitian matrix integral is proportional to e^ ^ and it thereby produces a solution which is 
singular at g = 0. More precisely, the loop equation l|4.25(l is very similar to that of the Kontsevich- 
Penner model |69j-j71| whose solution is given by the singular branch. While their equations of motion 
are the same, the complex and Hermitian matrix models differ in the choice of boundary conditions. 
The situation here is in marked contrast to that in the absence of the external field, whereby the dif- 
ference between partition functions calculated in the large N limit with different integration domains is 
exponentially small and hence does not affect the solution at leading order N = 00. 
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It follows that the function W{z) is given everywhere in the complex 2;-plane by the 
contour integral 



. / ^ ./ '-^■;;-'^' . (D,i3) 

J 477 1 w — z y [w — bi)[w — 02) 

[biM 

The discontinuity equation ()D.12|) determines the solution flD.13|) uniquely up to terms 
which are regular at z = 0, and the large z behaviour ()D.9|) implies that the regular terms 
vanish. Substituting ()D.1|) into ()D.13|) and integrating along the cut [&i,&2] then yields 

w(z) = ±_VWEMSEM 

2 J z-X ^(A-6i)(A-62) 

ai 

The large z behaviour ()D.9|) also generates two boundary conditions which unambiguously 
determine the branch points bi and 62 of the function W{z). Expanding the right-hand 
side of flD.14|) for z ^ 00 we encounter a constant term, which must vanish, and a term 
proportional to ^, with known residue —1. This supplements ()D.14|) with the respective 
constraints 

61 + 62 = -2~g I dA ^^^^ 



v/(A-6i)(A-62) 

3 (6? + 6^) +2^2 + 8^ = -8^ /dA^=M^L==. (D.15) 

J a/ A -61 A -62 



11 

12 



11 



Appendix E Calculation of the Convolution Kernels 

The trace formula ()5.6p for the integration kernel Ql can be written in terms of the 
generating function (j3.8j) for the Landau wavefunctions as 

X 'Pui,u2{xl)'Pu^^u3{xl-i) ■ ■■'Puz,uA^i) ■ (E.l) 

This representation can be proven by substituting the definitions ()3.7p into ()E.1|) and 
repeatedly applying the integral identity 

^ e-\^\'' = i\ (E.2) 

1 TT 

to reduce (jE.lj) to ()5.(jj) . Inserting the explicit forms of the generating functions ()3.8j) into 
(jE.ljl then leaves a set of L coupled Gaussian integrals. The result of these integrations is 
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most efficiently presented by introducing an L dimensional vector notation involving the 
L X L shift matrix 





/ n 
f 


U . . . 


n 1 \ 
\ 








1 


... 









Tl = 





1 . . . 







(E.3) 




Vo 


... 


1 o) 












/ 


C \ 




= (6 


, . . . 




H 




(E.4) 










a/ 




e = 






= VBz 




(E.5) 



and the complex L-vectors 
where 



are the dimensionless rescalings of the complex coordinates (j3.1|) on by the magnetic 
length. 

After some careful inspection, it is easy to see that the integrated form can then be 
written as 

L-l 



2716 



detfl + r,)-^ e^'^^i'^+^^r'i^W-i . 



(E.6) 



This same result could also have been obtained directly from ()5.3j) by using the Fourier 
representation of the delta-function. The expression ()E.6|) as it stands is formal, because 
while for L odd all eigenvalues of 1 + are non-zero, for L even there is exactly one zero 
eigenvalue and the matrix 1 + F^ is singular. We will therefore regulate the expression 
()E.6|1 by replacing the shift matrix F^ with aF^ and at the end take the limit a — > 1. 
Thus we will instead compute 

L-l 



.XL 



1 



2n9 



hm det(l + aFi)-i ^^P(-^^^TT^ 



(E.7) 



The determinant in ()E.7|1 can be computed by expanding in minors along the first row 
to produce two triangular determinants and hence get 





1 


. 


. 


a 




a 


1 . 


. 





det(l + aFL) = 





a . 


. 










. 


. a 


1 



-a] 



(E.8) 



The calculation of the inverse matrix in (jE.7jl is also straightforward owing to the identity 
1, which yields 



oo ^ L— 1 



(E.9) 
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Since the p^^ power of the shift matrix generates a cychc permutation of order p when 
acting on L-vectors, i.e. 

(riO/ = W)^... , (E.io) 

after a reshuffling of indices the quadratic form in ()E.7|) can be written as 



+ 2 5^(-a)^-^(0O + (-«)^O0 



KJ 



(E.ll) 



To get the final result for the kernel ()E.7|) . the cases of L odd and of L even should now 
be treated separately. 

When L = 2r + 1 is odd, the limit a — 1 is non-singular, and hence we can simply 
set a = 1 in ()E.8|) and ()E.11|) to compute ()E.7|) as 

g2r+i{xi,...,X2r+i) = (j^ exp ^- i ^,{-1)' X i ■ Bx^ . (E.12) 

When L = 2r is even, the regulator a must be removed carefully. For this, we set a = 1—e, 
and expand ()E.8|) and ()E.11|) to second order in e | to get 

det(l + ar2r) = 2re + 0{e^) , (E.13) 



2i 

r 



^(-1)^"^ {J-I~r)xi- Bxj + 0{e) . (E.14) 



KJ 



The first term in ()E.14jl combines with ()E.13jl to produce a delta-function in the limit 
e I owing to the identity 

lim — e~*l"l' = 6^^\x) . (E.15) 
elO 27ie 

This leads to the expression ()5.7|) quoted in the main text. 



Appendix F Calculation of the Free Propagator 

In this appendix we shall compute the propagator ()5.1ip . We will start by evaluating the 
Landau heat kernel 

oo 

Pt{x,y) ^ {x\e-'^'\y) = ^ e'^^^^"^)* M^) Mv) (F-l) 

e,m=i 
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for < t < cxD. For this, we will find it convenient to compute the function 

Xt{k,x) = e*°' e^'^-^ 0o(x) = J d'y Pt{x,y) e''-y Mv) , (F.2) 
where 0o is the ground state wavefunction ()3.6|) . By a simple calculation we find 

Xt=e'^' e^(^+^^)0o (F.3) 

where 

c = K3 + Kh, K = (^ki + i , (F-4) 

with a and b the ladder operators ()3.2|) . By using the Baker- Campbell-Hausdorff formula 
with [c, c''"] = 2|fi;p and c0o = we obtain 

=exp(i e*°' ct e-*°') e"!''!' e*"^' 0o • (F.5) 

From the oscillator representation (j3.16j) of the Landau Hamiltonian at B9 = 1 and (j3.22jl 

we get 

X, = e-2^* e-l'^l' exp i {k e"^^* aU^bt) 00 . (F.6) 

Inserting 0o = exp i (k e~^^* b + Ka)0o into ()F.6|) and using the Baker-Campbell- 
Hausdorff formula again yields 

XtiKx) = e-2^* e-l'^l' el'^l'^^P^-"^*) exp i [k e"^^* l + ^i) M^) , (F.7) 
where we have used (|3.2|) and the rescaling ()E.5|) . 



Along with ()F.2|) . the expression ()F.7|) allows us to compute the heat kernel ()F.1|) via 
the Fourier transform 

Ptix^y) = TTT e-l'^l^(i-P(-^^*)) 
0o(2/) i (2vr)^ 

X exp i [k e-^^* ^ + 7?^) . (F.8) 

The integral in ()F.8|) is Gaussian, and after a straightforward computation we arrive at 
the final result 

pj^^^y) = ^ ^ g-f coth(2Bt)|.-,P -i.-By /p_9^ 

^ 27rsinh(25t) ^ ^ 

Note that from ()E.15|1 it follows that Pt{x, y) 5^'^\x — y) for 5 | 0, as it should. Finally, 
the propagator ()5.1H) is given by 

oo 

C,{x,y) = jdt e-'^' {x\e-'^'\y) 



= — e-'^-^J' / dt e"*'^ ^ . (F.IO) 

27r J smh{2Bt) ^ ' 





The change of variables 



2 

then brings the propagator ()F.10|) into the form (|5.12|) given in the main text. 



u = ^ — ^ ( coth(25t) - 1 ) (F.ll) 
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Appendix G Scaling Limit of the Function 7^ (x, y) 



To compute the sum over Landau eigenf unctions ()5.18j) . it is convenient to introduce the 
generating function 

Proceeding as with the representation ()E.1|) . by using (|5.18p we may then express ()G.H1 
in terms of the generating functions ()H.7j) as 

^{x,y-t)=47r9 [ ^ e-H' / ^ P^,^ e--(x) e • (G.2) 

J ITC J ITI 



Upon substituting in the exphcit form (j3.8p . we observe that the integral over u is Gaus- 
sian, and that the resulting r integration yields the Bessel function Jq of the first kind of 
order 0, giving 



■t) = A ^-\--y?/'2e+i.-By+t^ J,(2t\x-y\/Vd^ . (G.3) 



By substituting the power series expansions of the exponential and Bessel functions into 
(lG.3j) and comparing with (jG.lll we find the functions ()5.18|) for finite C and 6 in the form 

.,(.,„^4,.-l,: e-l.-.l"/.»..... 'tjjf^. (^)' ■ 

We are interested in the function ()G.4|) in the scaling limit ^ oo, ^ — > oo with | 

fixed. In this limit, the exponential factor in ()G.3|1 reduces to e*^. The extraction of large 
orders in I in the Taylor series ()G.1|1 is accomplished by the standard method of contour 
integration to write 

/M 
^— . (G.5) 

i=0 

By rescaling the integrand of ()G.5|) as t = \/(. — \ w and using ()G.3|) we get 

dw 



7,(a;,y) = 4(£-l)! i jJ2w\x-y\^{l-l)/e) e (^"^H'^'-^i"-) . (G.6) 



w=0 



In the large I limit, the integral (jG.6j) can be evaluated by using the saddle-point approx- 
imation. There are two saddle points, at = ± 1. By expanding ()G.6|) near these saddle 



points and using the Stirling approximation l\ y2Txi e ^ £^ for ^ — > oo, we finally find 

7,(x,2/)=4Jo(2|x-y| , (G.7) 

which on trading the integers £ for the eigenvalues \i defined in fl4.31|l yields the result 
()5.19|1 in the main text. 



65 



References 



[1] A. Konechny and A. Schwarz, "Introduction to M(atrix) theory and noncommutative 
geometry," Phys. Kept. 360, 353 (2002) |arXiv:hep-th/0012145| , 0107251]. 

[2] M. R. Douglas and N. A. Nekrasov, "N oncommutative field theory," Rev. Mod. Phys. 
73, 977 (2002) |arXiv:hep-th7o T06048 1 . 

[3] R. J. Szab o, "Quantum field t heory on noncommutative spaces," Phys. Rept. 378, 
207 (2003) ^arXiv:hep-th70T09 162|. 

[4] S. Minwalla, M. Van Raamsdonk and N. Seiberg, "Noncommutative perturbative 
dynamics," J. High Energy Phys. 0002, 020 (2000) |arXiv:hep-th/9912072j. 

[5] H. O. Girotti, M. Gomes, V. O. Rivelles and A. J. da Silva, "The noncommuta- 
tive supersymmetric nonlinear sigma model," Int. J. Mod. Phys. A 17, 1503 (2002) 
|arXiv:hep-th/0102101|. 

[6] E. T. Akhmedov, P. De Boer and G. W. Semenoff, "Noncommutative Gross-Neveu 
model at large iV," J. High Energy Phys. 0106, 009 (2001) |arXiv:hep-th/0103199|. 

[7] I. Chepelev and R. Roiban, "Renormalization of quantum field theories on 
noncommutative R'^ I: Scalars," J. High Energy Phys. 0005, 037 (2000) 
(arXiv:hep-th/9 911098|. 

[8] I. Chepelev and R. Roiban, "Convergence theorem for noncommutative Feyn- 
man graphs and renormalization," J. High Energy Phys. 0103, 001 (2001) 
[arX iv:hep-th/000809 1 . 

[9] H. Grosse and R. Wulkenhaar, "Renormalization of $^-theory on noncommutative 
in the matrix base," arXiv:hep-th/0307017, 

[10] E. Langmann, R. J. Szabo and K. Zarembo, "Exact solution of noncommutative field 
theory in background magnetic fields," arXiv:hep-th/0303082 

[11] Y. Habara, "A new approach to scalar field theory on noncommutative space," Prog. 
Theor. Phys. 107, 211 (2002) arXiv:hep-th/0 105002 1. 

[12] E. Langmann and R. J. Szabo, "Duality in s calar field th eory on noncommutative 
phase spaces," Phys. Lett. B 533, 168 (2002) |arXiv:hep-th /0202039|. 

[13] E. Langmann, "Interacting fermions on noncommutative spaces: Exactly solvable 
quantum field theories in 2n + 1 dimensions," Nucl. Phys. B 654, 404 (2003) 
|arXiv:hep-th/0205287|. 

[14] Y.-G. Miao, H. J. W. Miiller-Kirsten and D.-K. Park, "Chiral bosons in noncommu- 
tative spacetime," arXiv:hep-th/0306034 

[15] C. Duval and P. A. Horvathy, "The 'Peierls substitution' and the exotic Galilei 
group," Phys. Lett. B 479, 284 (2000) arXiv:hep-th/0002233 . 

[16] V. P. Nair and A. P. Polychronakos, "Quantum mechanics on the noncommutative 
plane and sphere," Phys. Lett. B 505, 267 (2001) |arXiv:hep-th/0011172|. 

[17] B. Morariu and A. P. Polychronakos, " Quantum me chanics on the noncommutative 
torus," Nucl. Phys. B 610, 531 (2001) |arXiv:hep-th/0102157| . 



66 



J. Gamboa, M. Loewe, F. Mendez and J. C. Rojas, "The Landau problem 
an d no ncommutative quantum mechanics," Mod. Phys. Lett. A 16, 2075 (2001) 
(arXiv:hep-th/0104224|. 

A. Jellal, "Orbital magnetism of two-dimension noncommutative confined system," 
J. Phys. A 34, 10159 (2001) |arXiv:hep-th/0105303|. 



C. Duval and P. A. Horvathy, "Exotic Galilean symmetry in the noncom- 
mutative plane and the quantum Hall effect," J. Phys. A 34, 10097 (2001) 
[ar Xiv:hep-th/0106089j . 

S. Bellucci, A. Nersessian and C. Sochichiu, "Two phases of the noncommutative 
quantum mechanics," Phys. Lett. B 522, 345 (2001) |ar Xiv:hep-th/010 6138|. 

P. A. Horvathy, "The noncommutative Landau problem and the Peierls substitution," 
Ann. Phys. 299, 128 (2002) |arXiv:hep-th/0201007|. 

D. Bak, S. K. Kim, K.-S. Soh and J. H. Yee, "Exact wavefunctions in a noncommu- 
tative field theory," Phys. Rev. Lett. 85, 3087 (2000) arXiv:hep-th/0005253|. 

D. Bak, S. K. Kim, K.-S. Soh and J. H. Yee, "Noncommutative field the- 
ories and smooth^ commutative limits," Phys. Rev. D 63, 047701 (2000) 

|arXiv:hep-th/0 006087l . 

J.-M. Gracia-Bondi'a and J. C. Varilly, "Algebras of distributions suitable for phase- 
space quantum mechanics I," J. Math. Phys. 29, 869 (1988). 

G. Mandal, S.-J. Rey and S. R. Wadia, "Quantum aspects of GMS solutions of 
noncommutative field theory and large N limit of matrix models," Eur. Phys. J. C 
24, 495 (2002) |arXiv:hep-th/0111059|. 

C. Becchi, S. Giusto and C. Imbimbo, "The Wilson-Polchinski renormalization group 
equation in the planar limit," Nucl. Phys. B 633, 250 (2002) |arXiv:hep-th7 0202155j . 

F. Lizzi, P. Vitale and A. Zampini, "The fuzzy disc," ^arXiv:hep-th / 03 062471 

J. Ambj0rn, Y. M. Makeenko, J. Nishimura and R. J. Szabo, "Finite matrix 
models of noncommutative gauge theory," J. High Energy Phys. 9911, 029 (1999) 
jarXiv:hep- th/9911041|. 

J. Ambj0rn, Y. M. Makeenko, J. Nishimura and R. J. Szabo, "Lattice gauge fields 
and discrete noncommutative Yang-Mills theory," J. High Energy Phys. 0005, 023 
(2000) arXiv:hep- th/0004147 . 

S. Kim, C.-k. Lee and K.-M. Lee, "Quantum field dynamics in a uniform magnetic 
field: Description using fields in oblique phase space," Phys. Rev. D 65, 045009 
(2002) ar Xiv:hep-t h/0110249 . 

N. Seiberg and E. Witten, "String theory and noncommutative geometry," J. High 
Energy Phys. 9909, 032 (1999) |arXiv:hep-th/9908142|. 

N. Seiberg, "A note on background independence in noncommutative gauge theories, 
matrix model and tachyon condensation," J. High Energy Phys. 0009, 003 (2000) 
iarXiv:hep-th/000801 3|. 

J. Ambj0rn and S. Catterall, "Stripes from (noncommutative) stars," Phys. Lett. B 
549, 253 (2002) |arXiv:hep-lat/0209106 . 



67 



W. Bietenholz, F. Hofheinz and J. Nishimura, "Noncommutative field theories be- 
yond perturbation theory," ,arXiv:hep-th/021225 8 

F. Lizzi, R. J. Szabo and A. Zampini, "Geometry of the gauge algebra in 
noncommutative Yang-Mills theory," J. High Energy Phys. 0108, 032 (2001) 
farXiv:hep-th/0107115| . 

G. 't Hooft, "A planar diagram theory for strong interactions," Nucl. Phys. B 72, 
461 (1974). 

E. Brezin, C. Itzykson, G. Parisi and J.-B. Zuber, "Planar diagrams," Commun. 
Math. Phys. 59, 35 (1978). 

G. Landi, F. Lizzi and R. J. Szabo, "From large matrices to the noncommutative 
torus," Commun. Math. Phys. 217, 181 (2001) |arXi v:hep-th/9912130j . 

R. Gopakumar, S. Minwalla and A. Strominger, "Noncommutative solitons," J. High 
Energy Phys. 0005, 020 (2000) |a rXiv:hep-th/0003160j . 

R. Gopakumar, M. Headrick and M. Spradlin, "On noncommutative multisolitons," 
Commun. Math. Phys. 233, 355 (2003) | arXiv:hep-th/0103256| . 

Harish-Chandra, "Differential operators on a semisimple Lie algebra," Am. J. Math. 
79, 87 (1957). 

C. Itzykson and J.-B. Zuber, "The planar approximation II," J. Math. Phys. 21, 411 
(1980). 

M. Adler and P. van Moerbeke, "String orthogonal polynomials, string equa- 
tions and two-Toda symmetries," Commun. Pure Appl. Math. 50, 241 (1997) 
iarXiv:hep- th/9706182|. 

K. Takasaki and T. Takebe, "Integrable hierarchies and dispersionless limit," Rev. 
Math. Phys. 7, 743 (1995) ,arXiv:hep-th/ 9405096,. 

S. Kharchev, A. Marshakov, A. Mironov, A. Morozov and A. Zabrodin, "Towards 
unified theory of 2D gravity," Nucl. Phys. B 380, 181 (1992) arXiv:hep-th/9201013. . 

M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (National 
Bureau of Standards, 1964). 

J. Ambj0rn, C. F. Kristjansen and Y. M. Makeenko, "Higher genus correlators for the 
complex matrix model," Mod. Phys. Lett. A 7, 3187 (1992) arXiv:hep-th/920702qj. 

J.-M. Drouffe and J.-B. Zuber, "Strong coupling and mean field methods in lattice 
gauge theories," Phys. Rept. 102, 1 (1983). 

W. Bietenholz, F. Hofheinz and J. Nishimura, "The renormalizability of gauge 
theory on a noncommutative plane," J. High Energy Phys. 0209, 009 (2002) 
^arXiv:hep-th/0203151j . 

J. C. Breckenridge, V. Elias and T. G. Steele, "Massless scalar field theory in a 
quantized spacetime," Class. Quant. Grav. 12, 637 (1995) ,arXiv: hep-th /9501108| . 

A. Matusis, L. Susskind and N. Toumbas, "The IR/UV connection in noncommuta- 
tive gauge theories," J. High Energy Phys. 0012, 002 (2000) arXiv:hep-th/0002075. . 



68 



V. V. Khoze and A. Travaglini, "Wilsonian effective actions and the IR/UV mix- 
ing in noncommutative gauge theories," J. High Energy Phys. 0101, 026 (2001) 
[arXiv:hep-th/0011218|. 

Y. M. Makeenko and I. Che pelev , "Supersymmetric matrix models and the meander 



problem," |ajXw:hep-th/9601139| 



G. W. Semenoff and R. J. Szabo, "Ferm ionic matrix models," Int. J. Mod. Phys. A 
12, 2135 (1997) |arXiv:hep -th/9605140|. 

J. Ambj0rn, Y. M. Makeenko and K. Zarembo, "Supersymmetric matrix models and 
branched polymers," Nucl. Phys. B 482, 660 (1996) arXiv:cond-mat / 9 606041| . 

C. -S. Chu and F. Zamora, "Manifest supersymmetry in noncommutative geometry," 
J. High Energy Phys. 0002, 022 (2000) arXi v:hep-th/9912153| . 

D. Klemm, S. Penati and L. Tamassia, "Non(anti)commutative superspace," Class. 
Quant. Grav. 20, 2905 (2003) [arXiv:hep-th/0104190j . 

R. Abbaspur, "Noncommutative supersymmetry in two dimensions," Int. J. Mod. 
Phys. A 18, 855 (2003) arXiv:hep-th/011 0005 . 

J. de Boer, P. A. Grassi and P. van Nieuwenhuizen, "Noncommutative superspace 
from string theory," arXiv:hep-th/0302078, 

H. Kawai, T. Kuroki and T. Morita, "Dijkgraaf-Vafa theory as large N reduction," 
|ar Xiv:hep-t h/0303210 

I. Chepelev and C. Ciocarlie, "A path integral approach to noncommu tative super- 
space," J. High Energy Phys. 0306, 031 (2003) |arXiv:hep-th/0304118| . 



N. Seiberg, "Noncommutative superspace, M = \ supersymmetry, field theory and 
string theory," J. High Energy Phys. 0306, 010 (2003) |arXi v:hep-th/0 305248|. 

R. Britto, B. Feng and S.-J. Rey, "Deformed superspace, M = \ supersymmetry and 
(non)renormalization theorems," arXiv:hep-th / 0306 2T5[ 



M. Hatsuda, S. Iso and H. U mets u, "No ncom mutative superspace, supermatrix and 
lowest Landau level," |arXiv:hep-th /0306251| 

S. Ferra ra, M. A. Lledo and . Macia, "Supersymmetry in noncommutative super- 
spaces,' 



arXiv:hep- th/0307039[ 



E. Langmann, "Ex actly solvable models for 2D correlated fermions," 
|arXiv:cond-mat /0206045| 

A. Connes and J. Lott, "Particle models and noncommutative geometry," Nucl. Phys. 
Proc. Suppl. 18B, 29 (1991). 

L. Chekhov and Y. M. Makeenko, "The multicritical Kontsevich-Penner model," 
Mod. Phys. Lett. A 7, 1223 (1992) iarXiv:hep-th/9201033j . 

L. Chekhov and Y. M. Makeenko, "A hint on the external field problem for matrix 
models," Phys. Lett. B 278, 271 (1992) arXiv:hep-th/9202006 . 

K. Zarembo and L. Chekhov, "Multicut solutions of the matrix Kontsevich-Penner 
model," Theor. Math. Phys. 93, 1328 (1992). 



69 



